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Abstract 

We consider nonlinear recurrences generated from cluster mutations applied to quivers that have the 
property of being cluster mutation-periodic with period 1. Such quivers were completely classified by 
Fordy and Marsh, who characterised them in terms of the skew-symmetric matrix that defines the quiver. 
The associated nonlinear recurrences are equivalent to birational maps, and we explain how these maps 
can be endowed with an invariant Poisson bracket and/or presymplectic structure. 

Upon applying the algebraic entropy test, we are led to a series of conjectures which imply that 
the entropy of the cluster maps can be determined from their tropical analogues, which leads to a sharp 
classification result. Only four special families of these maps should have zero entropy. These families are 
examined in detail, with many explicit examples given, and we show how they lead to discrete dynamics 
that is integrable in the Liouville- Arnold sense. 

Keywords: Integrable maps, Poisson algebra, cluster algebra, algebraic entropy, tropical, monodromy. 

1 Introduction 

Cluster algebras were first developed by Fomin and Zelevinsky more than a decade ago [TU] . Their structure 
arises in diverse parts of mathematics and theoretical physics, including Lie theory, quantum algebras, 
Teichmuller theory, discrete integrable systems and T- and Y-systems. One of the original motivations for 
cluster algebras came from a series of observations made by Michael Somos and others (see |16j). concerning 
nonlinear recurrence relations of the form 

Xn+N X n = F(x n+ i, . . . ,X n+ N-l), (1-1) 

where F is a polynomial in N — 1 variables. The original observation of Somos was that certain choices of 
F lead to integer sequences when all N initial values are chosen to be 1. This was explained by the further 
observation that for such special F the recurrence (jl.ip exhibits the Laurent phenomenon, meaning that all 
iterates are Laurent polynomials in the initial data with integer coefficients. One of the most well-known 
examples is the Somos-4 recurrence given by 

2-n+4 2-n = 2-n+3 a 'n+l "I" 2 'n+2i -2) 

which generates the sequence 1, 1, 1, 1, 2, 3, 7, 23, 59, 314, 1529, 8209, 83313, . . Q starting from four initial Is, 
while if the initial data x\, X2, £3, £4 are viewed as variables then the iterates x n belong to the Laurent 
polynomial ring Zfof 1 , x^ 1 , xf 1 , xf 1 ]. 

In this paper we consider recurrences of the general form 

Xn+NX n = Y[ *n+j + II X n+P ( L3 ) 
a j > a j < 
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where the indices in each product lie in the range 1 < j < N— 1, with the exponents («i ajv— l) forming an 

integer (N — l)-tuple which is palindromic, so that Oj = a^—j. The main purpose of this paper is to identify 
which recurrences of the form (jl.3l) can be regarded as finite-dimensional discrete integrable systems, in the 
sense of the standard Liouville- Arnold definition of integrability for maps [351 [37] ■ The latter requires that 
a map should preserve a Poisson bracket, as well as having sufficiently many first integrals that commute 
with respect to this bracket. 

A quiver is a graph consisting of a number of nodes together with arrows between the nodes. In [13j . 
Fordy and Marsh showed how recurrences of the form (|1.3I) arc produced from sequences of mutations in 
cluster algebras defined by quivers with a special periodicity property with respect to mutations. They 
define a quiver Q with N nodes to be cluster mutation-periodic with period m if fi m •• - /U2 • MiQ = p m Q, 
where denotes quiver mutation at node j and p denotes a cyclic permutation of the nodes. Associated 
with the quiver mutation there is a corresponding cluster mutation acting on a cluster x = (x\, . . . , Xjv) in 
a (coefficient-free) cluster algebra. In the period 1 case, m = 1, the action of a suitable ordered sequence of 
cluster mutations on cluster variables is precisely equivalent to iteration of a recurrence of the form (|1.3p . 
A complete classification of period 1 quivers is given in jT5]: any such quiver produces a recurrence, and 
conversely any recurrence of the form (|1.3[) determines a cluster mutation-periodic quiver with period 1. 

1.1 Outline of the paper 

In the next section, we briefly review how the recurrences (|1.3p , or the equivalent birational maps in dimension 
N , come about from cluster mutations. The main object is the N x N skew-symmetric integer matrix B 
corresponding to the quiver, which not only defines the exponents appearing in (|1.3[) . but also produces 
a presymplectic form which is invariant under the map; this is the two-form introduced in |18) . When 
dot B 7^ 0, the form is symplectic, so the map automatically has a nondegenerate Poisson bracket. The main 
result of section 2 is Theorem 12.61 which states that (even if det B = 0) it is always possible to reduce (|1.3[) 
to a symplectic map, possibly on a space of lower dimension. This provides us with the appropriate setting 
in which to consider Liouville integrability in the rest of the paper. 

In section [3] we consider the recurrences (|1.3j) in the light of the algebraic entropy test [2] . We give 
details of a series of conjectures which show that the algebraic entropy can be determined explicitly from 
the tropical version of (jl.3[) , expressed in terms of the max-plus algebra. From the point of view of the rest 
of the paper, the main result is a corollary of these conjectures (Theorem I3.12p . which classifies the cases 
with zero entropy into four families, labelled (i)-(iv). All of the maps in family (i) have periodic orbits, and 
it is a trivial task to show that they are Liouville integrable maps. The majority of the rest of the paper is 
devoted to families (ii), (iii) and (iv). 

Section U is concerned with the family (ii) , which arises from cluster mutations applied to the primitive 
quivers, denoted Pfi' , which are defined for each positive integer N and q = 1, . . . , L^V/ 2J . These were intro- 
duced in |13j , where they were shown to be the building blocks of all mutation-periodic quivers with period 
1. They are also equivalent to affine A- type Dynkin quivers (or copies of such), and it was shown by Fordy 
and Marsh that the cluster variables in this case satisfy linear recurrence relations with constant coefficients. 
(This was subsequently shown for the general case of cluster algebras associated with affine Dynkin quivers, 
in pQ and |23|.) Here we give a new proof of these linear recurrences, which relies on additional linear rela- 
tions with periodic coefficients, and their associated monodromy matrices. These periodic quantities are the 
key to the Liouville integrability of the maps in family (ii) . A large number of new examples of integrable 
maps arise in this construction and are explicitly presented. 

Section [5] deals with the family (iii), each member of which arises from a quiver which is a deformation of 
a primitive Pjf , for a particular q and N. The general properties of this family are very close to those of the 
primitives. In particular, the cluster variables satisfy linear recurrences with constant coefficients, and there 
are additional linear relations with periodic coefficients. Once again, associated monodromy arguments, and 
Poisson subalgebras defined by the periodic quantities, are the key to the Liouville integrability of members 
of this family. Again, many new examples of integrable maps arise and are explicitly presented. 

The family (iv) consists of Somos-type recurrences with three terms, typified by (|1.2j) . In section [6] we 
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outline some different approaches to understanding the Liouville integrability of the maps in this family, such 
as making reductions of the Hirota-Miwa equation and its Lax pair and by deriving higher bilinear relations 
with constant coefficients. 

Some of our results were announced previously in |15j . 

2 Symplectic maps from cluster recurrences 

Given a recurrence, a major problem is to find an appropriate symplectic or Poisson structure which is 
invariant under the action of the corresponding finite-dimensional map. Remarkably, in the case of the 
cluster recurrences f| 1 . 3[) this problem can be solved algorithmically. 

2.1 Recurrences from periodic quivers 

A quiver Q with TV nodes and no f-cycles or 2-cycles admits quiver mutation. The mutation at node k 
produces a new quiver Q = pkQ which is obtained as follows: (i) reverse all arrows in/out of node k; (ii) if 
there are p arrows from node j to node k, and q arrows from node k to node £, then add pq arrows from 
node j to node £; (iii) remove any 2-cycles created in step (ii). 

An TV x TV skew-symmetric matrix B with integer matrix elements bjg defines a quiver Q with TV nodes, 
without 1-cycles or 2-cycles. Matrix mutation applied at the vertex k is also denoted pk, and starting from 
B it produces a new matrix B — pkB = (bji) defined by 

- _ f -bjt if j = k or £ = k, ,^ ^ 

il \ b ]i + \{\b]k\bki + bj k \bki\) otherwise. 

This matrix mutation is equivalent to the action of fik on Q, via quiver mutation, producing the new quiver 
Q = PkQ- In addition to the matrix mutation, the cluster variables x = (xi, X2, ■ ■ ■ , xn) are transformed by 
Pk to a new cluster x = (ii, X2, ■ ■ ■ , xn) in such a way that all variables except Xk are left unaltered, so that 
xj = Xj, j =/= k, and the exchange relation corresponding to the mutation pk is conveniently written as 

N N 

■n...r, • • \[r '■ . (2.2) 

j=l j=l 

where [b] + = max(&, 0). Note the identity ^(o|&| — \a\b) — a[— b] + — b[—a] + . 

In what follows, we require that the N x N skew-symmetric matrix B defines an TV-node quiver Q that 
is cluster mutation-periodic with period 1. All such matrices were classified in |13j . Cluster mutation- 
periodicity, in the case that the period is 1, means that after applying a single step of mutation at one of the 
nodes, /zi say, the quiver Q is the same as the quiver pQ obtained from Q by applying the cyclic permutation 
p, given by p : (1,2,3,..., N) H> (TV, 1, 2, . . . , TV — 1), such that the number of arrows from j to k in Q is 
the same as the number of arrows from to p _1 (fc) in pQ. This periodicity requirement corresponds 

to explicit conditions on the matrix elements of B, namely that 

bj+i.k+i = bjk (2.3) 

for all j, k, where the indices are read modulo TV. Since p N = id, from the periodicity it is clear that Q is 
preserved by the composition of TV mutations that cycle around its nodes, i.e. p^ • ■ ■ fJ>2 ■ p\Q = Q- 
Starting from B, one constructs the TVth order recurrence relation 

2V-1 N-l 

x n+N x n = I] x[^f ll+ + x^ +l] + , n = 1, 2, 3, . . . . (2.4) 

3=1 3=1 

If B satisfies the conditions (|2.3p . then iterating the recurrence (|2.4[) starting from the initial data (xi, . . . , xn) 
is precisely equivalent to applying cluster mutation p\ at the vertex 1, followed by the subsequent mutations 
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/X2,A t 3i • • • iPjiMj+i) • ■ •) an d so on. The recurrence relation (|2.4[) is clearly reversible, in the sense that (as 
long as neither is zero) it can be solved both for x n +N, to iterate forwards, and for x„, to iterate backwards. 
This means that the index n in (|2.4I) is allowed to take all values n € Z, and also that iteration of the 
recurrence is equivalent to iteration of the birational map tp from to itself, defined by 



Xl 

x 2 



\ 



XJV-1 

V x N J 



( x 2 \ 

X3 



XN 
\ X N+ i J 



N-l Jb ljj+1 ]+ 



where 



XN+l 



(2.5) 



One can decompose the map (|2.5|) as ip = p^ 1 ■ p,i, where in terms of the cluster x = (xj) the map \i\ 
sends (xi, X2, ■ • . , x/v) to (xi, £2, • • • , xjv), with Xi defined according to the exchange relation (|2.2[) with 
k = 1, and sends (xi, X2, . . . , xn) to (X2, • • • , Xn-> x i)- Due to the periodicity requirement on B, we have 



P 



pi (B) = B, so the action of ip on this matrix is trivial. 



2.2 The Gekhtman-Shapiro-Vainshtein bracket 

In [T7] it was shown that very general cluster algebras admit a linear space of Poisson brackets of log-canonical 
type, compatible with the cluster maps generated by mutations, and having the form 



{Xj Xfc} — Cjk XjX]^ 



(2.6) 



for some skew-symmetric constant coefficient matrix C = (cjk)- Compatibility of the Poisson structure 
means that the cluster transformations pi given by (|2.2p correspond to a change of coordinates, x = /Xi(x), 
with their bracket also being log-canonical, 



{.Xj, Xfc} — Cjk XjXfc 7 



(2.7) 



for another skew-symmetric constant matrix C = (cjk)- 

Our viewpoint is to regard the cluster transformation as a birational map x 1— > x = ip (x) in C , and 
require a Poisson structure that is invariant with respect to ip (not just covariant). Therefore, in (|2.7[) we 
require C — C . However, there may not be a non-trivial log-canonical Poisson bracket that is covariant or 
invariant under cluster transformations. 



Example 2.1. Corresponding to (|2.5p with N — 3, the matrix B and birational map on C 3 are given by 

(2-f 





x 2 

(X2X3 + l)/xi 



B = 



Suppose that there is an invariant Poisson bracket of the form (|2.6|) . The condition <p*xj = Xj+i, implies 
that Cj+i.k+e — Cjk for all indices j, k,£ in the appropriate range, so C is a Toeplitz matrix: 



C = 



Upon taking the bracket of both sides of the relation x±x\ = X2X3 + I with x 2l one finds that {a — fi)x\XiXt± = 
—olx\xz, so (a — /3)(x2X3 + x 2 ) — —olx\xz, which gives a — j3 = 0, so the bracket is trivial. 

Remark 2.2 (Non-invariance of symplectic leaves). Even in the case where the map ip does admit an 
invariant log-canonical Poisson bracket, it may be degenerate, and in that case <p need not preserve the 






a 


13 


—a 





a 


-P 


—a 






4 



symplectic leaves of the bracket. For instance, the Somos-4 recurrence (|1.2[) has the invariant Poisson 
bracket [21] 

{xj,x k } = (k-j)xjXk, (2.9) 

which has rank two, but the two independent Casimirs x\x^jx\, X2X^jx\ are not fixed by the action of ip (see 
Example 12.111 below), so the symplectic leaves are not preserved. The analogous observation for Somos-5 
appears in [15] . 

In general, we shall see that it is more useful to start with a two-form in the variables Xj, rather than a 
Poisson bivector field corresponding to a bracket. 

2.3 Symplectic forms for cluster maps 

Given a skew-symmetric matrix B, one can define the log-canonical two-form 

lu = J2~^- dx 3 ;Ada; fc , (2.10) 
j<k XjXk 

which is just the constant skew- form uj — Ylj< k bjk dzj A dz k , written in the logarithmic coordinates 
zj = log Xj , so it is evidently closed, but may be degenerate. In [TJ] (see also [5]) it was shown that for a cluster 
algebra defined by a skew-symmetric integer matrix B, this two- form is compatible with cluster transforma- 
tions, in the sense that under a mutation map \ii : x t— > x, it transforms as /i*uj = X^<fc bj k d\og£j Adlogifc. 

In the case that the matrix B is nondegenerate, w turns out to be a symplectic form for the map ip, but 
in general it is a presymplectic form. For the purposes of our discussion it is convenient to present some 
formulae from 1131 . which give the following result. 

Lemma 2.3. Let B be a skew- symmetric integer matrix. The following conditions are equivalent. 

(i) The matrix B defines a cluster mutation-periodic quiver with period 1. 

(ii) The matrix elements B satisfy the relations 

b JtN = bij+i, j = l,...,N-l, (2.11) 

and 

bj+i,k+i = b jtk + b lij+1 [-bx,k+i]+ - b 1 ,k+i[-b 1 j +1 } + , l<j,k<N-l. (2.12) 

(iii) The two-form w is preserved by the map ip, i.e. (p*ui = uj. 

Proof: The proof of (i) (ii) follows straight from the definition of periodicity (see the proof of Theorem 
6.1 in [IS])- The implication (ii) (iii) in Lemma [2.31 is a consequence of Theorem 2.1 in [18] , after 
applying the permutation p _1 to the coordinates. The reverse implication follows from a direct calculation, 
which is omitted. □ 

The formulae (12.11[) and (|2.12[) entail that for a period 1 cluster mutation-periodic quiver, the matrix 
B is completely determined by the elements in its first row, so that each recurrence of the form (|1.3p with 
palindromic exponents corresponds to a matrix B. Theorem 6.1 in [13j is equivalent to the following formula 
for b jk : 

max(j-l,A r -fe) 

bjk = -bkj = a k -j + ^ at[-ai+k-j]+ - at+k-j[-at]+, 1 < j < k < N, (2-13) 

t=i 

where the a,j = for j = 1, . . . , N — 1 form a palindromic integer (JV — l)-tuple, such that a,j — a^-j. 

Apart from being skew-symmetric, B is also symmetric about the skew diagonal, i.e. bjk = &jv-fc+i,iv— j+i- 
Henceforth when we refer to a recurrence of the form (|2.4p . and the corresponding matrix B — {bjk), we 
assume that the conditions (|2.11|) and (|2.12|) hold. 
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Example 2.4 (Corollary 2.2 in 21). For integer values of c > 0, the skew-symmetric matrix 

(2.14) 



B = 



( -1 c -l\ 

1 -(c+1) c 

-c (c+1) -1 

\ 1 -c 1 / 



satisfies the conditions (I2.11j) and (|2.12l) . which means that it defines a period 1 cluster mutation-periodic 
quiver Q, in the sense of |13) . Thus, by Lemma 12.31 for each c the map ip corresponding to the recurrence 

(2.15) 

preserves the two-form 

/ dxi A da;2 dxi A dx4 dx% A dx4 \ / dxi A d^3 da^ A dx^ \ dx2 A dx^ 

' -(c+1) . (2.16) 



\ XiX 2 X1X4 X 3 X4 ) \ X\X% X 2 X 4 J X2X3 

Remark 2.5 (Invariant volume form). For all maps ip of the form (|2.5p . or of the more general form 
the volume iV-form 

»= d " lA w -- AdXjV , (2.17) 

Uj = l X 3 

is invariant up to a sign, depending on the parity of N, i.e. (p*fl = ( — 1) 0. In the case that ui is 
nondegenerate, which can only happen for even N — 2K, up to overall scale this volume form is the 
Poincare invariant A ... A w (K terms). 

If the matrix B is nondegenerate, then (up to rescaling by an overall constant) the associated log-canonical 
Poisson bracket (|2.6|) is given by the dual bivector field 



J = J2 c 3k XjX k g^^Q^ with c = (cj*) = B-\ (2.18) 

j<k 



In the case that det B = 0, it is necessary to consider a projection to a lower-dimensional space with a 
symplectic form, as follows. 

Theorem 2.6. The map ip is symplectic whenever B is nondegenerate. For rank£> = 2K < N, there is a 
rational map it and a symplectic birational map (p such that the diagram 

C N — £-> C N 
1" 1" (2 ' 19) 

£ 2K C 2K 

is commutative, with a log-canonical symplectic form Cj on C 2K that satisfies 7r*tli = ui. 
The proof of this theorem will occupy most of the rest of this subsection. 

To begin with we consider the null distribution of oj, which (away from the hyperplanes Xj = 0) is 
generated by N — 2K independent commuting vector fields, each of which is of the form 

d N d 

= — for u = (uj) e ker_B. (2.20) 



dt ^ J J dx 

7=1 



Since this is an integrable distribution, Frobenius' theorem gives local coordinates t\, . . . , tN-2K, yi, ■ ■ ■ , yiK 
such that the integral manifolds of the null distribution are given by yj =constant, j — 1,...,2K. The 
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coordinates yj must be invariants for these commuting vector fields, and can be chosen as linear functions of 
the logarithmic coordinates Zj — log xj , but for our purposes it is more convenient to take functions of the 
form 

y = x v '-=Y[ x 7 < * => (u,v)=0 VuekerB, 

j 

where ( , ) denotes the standard scalar product. This yields a log-canonical symplectic form in terms of yj, 
which is the generalised Weil-Petersson form in [T5] . 

Alternatively, one can consider integral curves of the vector fields (I2.20[) . each of which is the orbit of the 
point x under the action of the algebraic torus C* , which acts by scaling the coordinates Xj . We denote this 
one-parameter group action by 

x4A".x=(A"'i j ), AeC*. (2.21) 

Combining N — 2K independent vector fields of the form (|2.20p , which (without loss of generality) can be 
defined by choosing N — 2K independent vectors u = (uj) € ker B with components Uj £ Z, we see that the 
integral manifold through a generic point x is the same as its orbit under the scaling action of the algebraic 
torus (C*)^ -2 ^. The coordinates yj = x Vj are the invariants under these scaling transformations, and, by 
choosing the Vj to be integer vectors, define a rational map it. 

Lemma 2.7. Suppose that the integer vectors Vi,...,V2k form a basis for imB, and define the rational 
map 

tt: — ► C 2K 

xi — > y, 2/j=x v ^, j = l,...,2K. 

Then the yj are a complete set of Laurent monomial invariants for all of the scaling transformations 112.21]) 
defined by integer vectors u G keri?, and there is a log-canonical symplectic form 

Cj = J ' fc dyj A dyk , with n*ui = uj. (2.22) 

Proof: Since the integer matrix B is skew-symmetric, the vector space has an orthogonal direct sum 
decomposition Q N =imB© ker B. The scaling action on Laurent monomials gives A u • x v = . \ UjVj Xj J — 
/ \( u > v ) x v , hence x v is invariant under the overall action of (C*) w - 2 ^ if anc l nly if (u, v) = for all u g ker B, 
so v G imB, and a basis of imB gives 2K independent monomial invariants. Now extend the basis of imB 
to a basis {vj}^ =1 for Q N , so that Vj £ kerf? for 2K +1 < j < N, and let M be the matrix whose rows 
consist of the N basis vectors. Then one finds the block matrix 

B^ = M~ T BM~ l = 

(with M~ T = (M _1 ) T , and T denoting transpose) where B = (bjk) is a nondegenerate skew-symmetric 
2K x 2K matrix. Defining the two-form (|2.22[) in terms of B gives tt*u} — u>, as required. □ 

The map tt and the specific form of w depend on the choice of integer basis for im_B. Here we consider 
rational maps between fixed affine spaces, even if these are not defined everywhere. In the case where 
B is nondegenerate, the diagram (|2 . 19[) is trivial for tt = id, but there are other non-trivial choices of tt, 
corresponding to different integer bases for imB. 

The following classical theorem (which is a special case of Theorem IV. 1 in [29]) provides a canonical 
choice of 7r, and via Lemma 12.71 gives Darboux coordinates for the presymplectic form uj. The proof in [29] 
is constructive. 

Theorem 2.8. If B is a skew- symmetric matrix of rank 2K in Mat at (Z) then there are integers hi, ... , hx 
and a unit matrix M € Matjy(Z) such that B = M T B^M , where 

B * = h i ( _°i J ) ® h * ( _°! J ) © • • • ® ( _°i I ) © °. 
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with hj\hj + i for j = 1, . . . , K — 1. 



In the above result, the first 2K rows of the unit matrix M provide a Z-basis for the Z-module imBz = 
imBfl 7L N (a sublattice of Z^), and the remaining N — 2K rows provide a Z-basis for keri?z. For our 
purposes, it is only the choice of integer vectors spanning im B that matters for the definition of the map n, 
and if the integers hj are not all 1 then it is not essential to have a Z-basis for va\B%. 

The second part of the proof of Theorem 12.61 involves showing that, when B is degenerate, it is possible 
to choose integer vectors v.,- spanning im_B in a way that is compatible with the map tp. In other words, 
ip should reduce to a symplectic map (p in the coordinates yj, and we require that the latter map should 
also be birational. The next result provides two sets of sufficient conditions on the vectors Vj which ensure 
birationality of the map (p. 

Lemma 2.9. Let {v,-}|ffi be integer vectors that span imB, and suppose that the columns of B belong to 
< Vi, . . . , v 2 isr >z- If either 

(a) imB z =< v i; . . . , v 2 k >z, or 

(b) each Vj belongs to the Z-span of the columns of B, 

then for 2/j(x) = x Vj , there exist rational functions fj(y), f](y) € Q(y) = Qd/ii ■ ■ • ; Vik) such that 
^*2/j( x ) = n*fj(y), (iy9~ 1 )*y J -(x) = 7r*/j(y) for j = 1,...,2K, and the map $ : y (fj(y)) is birational 
and symplectic. 

Proof: Let y(x) = x v be one of the coordinate functions yj in the image of the map tt, with r\ = logy, and 
denote the columns of B by Wi, . . . , Wjy. Then 



where m is an integer vector with components mi = — ujy, TUj = Vj-i + [— bij]+vrsr for 2 < j < N. 

By definition the Wj span imB, so we can write v = Ylk—i dk w k for some dk & Q, which in components 
gives Vj = X^fcLi dkbjk for j = 1, . . . , iV. By using (|2.1ip and (|2.12[) one obtains 




<^*f?(x) = ^w fe <,5*2: fc = ^ v k z k +i + v N (-Zi + y^[-&ij]+gj +log(l + expf) 

fc=l fc=l 3=2 



with i* 1 given by _F := X^fcL] 1 ^l.fe+i-Zfc+i- By the initial assumption on the Vj, we can write wi = Y^ij=i c j w j 
for Cj 6 Z, which implies that expi* 1 = tt* expj^ ■ CjT)j = tt* . , the pullback of a rational function of y. 
Then wc have 



tp*ri(x) = (m, z) + log(l + exp F) VN , (2.23) 



riij 




for j ' > 2, while mi 



Lf=i <4&JVfc = Y>k=i d kh,k+l, implying that 



JV-l 



N 




Hence we see that m £ \Yt\B%. In case (a), m 
exponentiating yields ip*y(x) — tt* f with 



Sj=i f° r ^7 € Z, so substituting in (12.23[) and 
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as required. In case (b), on the other hand, each dk € Z and the are in < Vi, . . . , v 2K >z, so again one 
can expand m in this basis with coefficients cj G Z, and the same formula holds for /. For the inverse map 
tp- 1 : (xi, . . .,x N ) n- (xo,X!, . . .,x N -i), with 



x 

X N 




in 



the fact that (ip 1 )*y(x) = n* p with /^(y) € Q(y) follows from an almost identical calculation. 

The rational map (p is defined by the 2K functions fj(y), with a rational inverse (p~ x given by the 
functions /J(y), and by construction ip ■ tt = tt ■ ip. Combining part (iii) of Lemma 12.31 together with the 
result of Lemma \2.7\ gives <p*ui — to = ip*Tt*Qj — tx*Cj = Tr*(ip>*u} — &) = 0, so (p*Q) = was required. Thus the 
proof of the lemma and the proof of Theorem 12.61 are complete. □ 



Remark 2.10. Theorem [2J3 provides a basis of im B which satisfies condition (a) above. To satisfy condition 
(b), one can choose any 2K independent columns (or rows) of B, corresponding to the r-coordinates of 

nansj. 



To illustrate Theorem 12.61 we now present a couple of examples. 

Example 2.11 (Somos-4). The Somos-4 recurrence (|1.2[) is the special case c = 2 of (|2.15|) . in which case 
the skew-symmetric matrix B of (|2. 14|) is degenerate, and ker_B is spanned by the integer column vectors 
Ui = (1,1,1,1) T , u 2 = (1,2,3,4) T . By (|2.21[) . each of these vectors generates a scaling action on the 
phase space C 4 , with weights given by their components, so that the torus (C*) 2 acts via (#1,2:2, £3, £4) — > 
(AiA 2 xi, A1A2 x 2l AiA| x 3 , AiA| Xi), for (Ai,A 2 ) e (C*) 2 . Then imB = (kcrB)- 1 =< vi,v 2 >, with vi = 
(1, —2, 1,0) T , v 2 = (0, 1, —2, 1) T , whose components provide the exponents for the monomial invariants 

2/1 = — 2/2 = — —■ (2.24) 

x 2 x 3 

(These monomials also provide two independent Casimirs for the degenerate Poisson bracket (|2.9p mentioned 
above. The coefficients of (|2 .9[) are obtained from the bivector Ui A u 2 .) This defines the rational map 
7r : C 4 — > C 2 from the Xj to the yj. Upon computing the pullback of ip on these monomials, one has the map 

* : (to+iW*)- (2 - 25 > 

which is of QRT type |33j , and preserves the symplectic form 

uj = dy 2 A dyi (2.26) 
2/12/2 

where u) — tt*lu, with u) given by the formula (j2 . 16[) for c = 2. 

Note that, in the preceding example, the chosen basis is such that every integer vector in imB can be 
written as a Z- linear combination of the vectors vi, V2, so the condition (a) in Lemma [2.91 holds. At the 
same time, from the matrix (|2.14[) with c — 2 we see that the first column of B is v 2 , and the last column is 
— vi, so condition (b) holds as well. The next example shows that (a) and (b) can give inequivalent results. 

Example 2.12. In [5T], a singularity confinement pattern was used to obtain the sixth-order recurrence 

X n +& X n — [x n +b x n+l) + x n +4 X n+3 X n+2 1 (2.27) 



9 



which is associated with mutations of a skew-symmetric matrix of rank 2, namely 



B 



( o 

2 

-2 
-4 
-2 
V 2 



-2 

G 
6 


-2 



2 4 2 

-6 -6 

-6 -6 

6 0-6 

6 6 

-4 -2 2 



-2 \ 
2 
4 
2 

-2 

o y 



A basis of imi? corresponding to case (a) of Lemma l2.9l is given in terms of the first and last columns of B 
by vi = — |wg, V2 = |wi, which gives the y-coordinates j/j = ^^^/(xj+ia^^rEj+a), j = 1, 2, and (up to 
rescaling) produces the same symplectic form as in (|2.26[) . The corresponding symplectic map is 



<P 



Hi 

m 



y-2 



(yi + i)/(waa) 



(2.28) 



W6, 
and 



whose singularity pattern under successive blowups was considered by Viallet [36 . The map (|2.28[) has 
positive algebraic entropy, indicating nonintegrability (See Example 13.71 in the next section.) 

However, one can take a different basis, corresponding to case (b) of Lemma 12.91 given by Vj = 
v 2 — Wi , which is not a Z-basis for im.B%. From this basis, one has the coordinates (2/1,2/2) = (2/1, 2/1 
the map becomes 

v ■ U2 y v (2/2 + 1)7(2^2) 

The map from (2/1, 2/2) to (y[, y' 2 ) is ramified, as generically there are four pairs of values (±2/1 > +2/2) for each 
pair (2/1,2/2)1 an( i so the two maps tp and are not conjugate to each other. 



3 Algebraic entropy and tropical recurrences 

The deep connection between the integrability of maps and various weak growth properties of the iterates 
has been appreciated for some time (see [37] and references). In the case of rational maps, Bellon and Viallet 
[2] considered the growth of degrees of iterates, and used this to define a notion of algebraic entropy. Each 
component of a rational map tp in affine space is a rational function of the coordinates, and the degree of 
the map, d — deg tp, is the maximum of the degrees of the components. By iterating the map n times one 
gets a sequence of rational functions whose degrees grow generically like d n . At the nth step one can set 
d n = deg tp n , and then the algebraic entropy £ of the map is defined to be £ = lim n _>. 0o i log d n . Generically, 
for a map of degree d, the entropy is log d > 0, but for special maps there can be cancellations in the rational 
functions that appear under iteration, which means that the entropy is smaller than expected. 

It is conjectured that Liouville- Arnold integrability corresponds to zero algebraic entropy. In an algebro- 
geometric setting, there are plausible arguments which indicate that zero entropy should be a necessary 
condition for integrability in the Liouville- Arnold sense [3]. In the latter setting, each iteration of the map 
corresponds to a translation on an Abelian variety (the level set of the first integrals), and the degree is a 
logarithmic height function, which necessarily grows like d n ~ Cn 2 . 

Often the algebraic entropy of a map can only be guessed experimentally, by calculating the degree 
sequence (d n ) up to some large n and doing a numerical fit to a linear recurrence. This is increasingly 
impractical as the dimension increases, and provides no proof that the linear relation, with its corresponding 
entropy value, is correct. In dimension two, exact results are possible via intersection theory [351 136) . 

In the rest of this section we seek to isolate those recurrences with £ = 0, by finding a condition on the 
exponents which should be necessary and sufficient for £ > 0. The main conjecture is the following. 
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Conjecture 3.1. For a birational map given by i2.5\) , corresponding to a recurrence of the form \2.J$ , the 
braic entropy is £ — log |A ma:E |, where, of all the roots of the two polynomials 

N-l 

P±(X) = X N + 1-J2 [±&i J+ i]+A J , (3.1) 
i=i 

is the one of largest magnitude. The entropy is positive if and only if 

(N-l N-l \ 

E l bl d+i]+ > E > 3 > ( 3 - 2 ) 

3=1 3=1 J 

We now give very strong evidence for the above assertion, showing how it rests on a sequence of other 
conjectures. For the recurrences (|2.4I) . the key to calculating their entropy is the Laurent phenomenon, which 
leads to an exact recursion for the degrees of the denominators. The Laurent property for the associated 
cluster algebra [11] implies that the iterates have the factorized form 

N ( ) ^ 
x n = " , v with N„ G Z[x] = Z[xt,...,x N ], M„ = x j " , 
-^t n (xj 

3=1 

where the polynomials N„ are not divisible by Xj for 1 < j < N, and M ra are Laurent monomials. A lower 

bound for the entropy is provided by the growth of degrees of denominators, and if the exponents dn are 
all positive (for large enough n) then the monomial M„ is the denominator of the Laurent polynomial x n , 
and N„ is the numerator. 

In addition to being Laurent polynomials, the form of the exchange relations (|2 .2[) in a cluster algebra 
means that all the cluster variables are given by subtraction-free rational expressions in terms of any initial 
cluster x = (x±, . . . , xn). This implies that the dynamics of the M n can be decoupled from the N„. 

Proposition 3.2. For all n, the exponent dn of each variable in the Laurent monomial M„ satisfies the 
recurrence 

(N-l N-l \ 

[b lyJ+1 ] + d n+j , [— 6i,j + i] + d n+ j , (3.3) 
3=1 3=1 / 

with the initial conditions d\ = —1, d 2 = . . . = e?jv = (up to shifting the index). 

Proof: Upon substituting x n = N„/M„ into (|2.4[) and comparing monomial factors on both sides one has 

(N-l N-l \ 

n M n+r il+ . n M L+r +il+ . ( 3 - 4 ) 
3=1 3=1 / 

where 1cm denotes the lowest common multiple. To be more precise, for any sequence of Laurent polynomials 
(x n ) with positive or negative coefficients, the formula (I3.4[) certainly holds provided that the two products 
on the right are not of equal degree in any of the variables Xj, for 1 < j < N. If it happens that D := 

- Y^t=i 1 [-°i,e+i]+ d n+£ for some h tnen tne coefficient of Xj D in each of the two terms 
on the right hand side of (|2 .4[) is a non-zero subtraction-free rational expression in the other variables Xk 
with l<k<N,k^j, and the sum of these two coefficients cannot vanish. Hence no cancellations can 
occur between the two products of numerators N„ on the right, and the formula (|3.4p always holds. Taking 
the degree of each variable on the left and right hand sides of (|3.4p gives the same recurrence (|3. 31) in each 
case. From the initial data x\, . . . , xn for (|2.4p it is clear that d^ = —1 and d^ = for n = 2, . . . , N, and 
the degrees dn for 2 < j < N have the same initial data but shifted along by an appropriate number of 
steps. □ 
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Remark 3.3. The recurrence p. 31) is the tropical (or ultradiscrete [3D]) analogue of the original nonlinear 
recurrence (|2.4I) . in terms of the max-plus algebra. It is a special case of the recursion for the denominator 
vectors in a general cluster algebra, which is stated without justification as equation (7.7) in [T2"] . 

Conjecture 3.4. Suppose that the sequence d n is not periodic and satisfies i3.S\) with the initial conditions 
as in Proposition^^ Then 

(a) d n > for all n > N , and 

(b) the total degree of the numerators satisfies degN„ ~ Gd n as n — > oo, for some constant C > 0. 

Remark 3.5. Part (a) above follows from the first part of Conjecture 7.4 in [T^]. Part (b) implies that the 
growth of the denominators completely determines the algebraic entropy, since the numerators grow at the 
same rate; it should be a consequence of Proposition 6.1 in [12] (graded homogeneity of cluster variables). 

Example 3.6 (Tropical Somos-4). The tropical version of the Somos-4 recurrence is 

d n+4: + dn = max(d n+3 + d n+1 , 2d n+2 ). (3.5) 

With initial conditions d\ = — 1 and d 2 = d 3 = d^ = this generates a sequence that begins 

-1, 0, 0, 0, 1, 1, 2, 3, 3, 5, 6, 7, 9, 10, 12, 14, 15, 18, 20, 22, 25, 27, 30, 33, ... , 

which are the degrees (in each of the variables xi,X2,xa,X4,) of the denominators of the Laurent polynomials 
generated by (|1.2[) . The preceding sequence has quadratic growth, d n ~ Cn 2 as n oo (consistent with 
the growth of logarithmic height on an elliptic curve), so that the algebraic entropy is zero, and this can be 
proved by considering the combination 

Y n = d n+ 2 — 2d n+ i + d n , (3-6) 

whose coefficients are the exponents in (]2.24[> . The sequence of quantities Y n is generated by the tropical 
analogue of the map (12.25[) . that is 

\Yi ) \ [Yi] + - 2Y 2 - Y 1 

and all of the orbits of the latter are periodic with period 8 (which is a special case of Nobe's results on 
tropical QRT maps [30]). Thus, if S denotes the shift operator such that Sf n = f n +i for any function of 
the index n, then applying the operator S 8 — 1 to both sides of (|3 . 6[) implies that the sequence of degrees 
d n satisfies a linear relation of order 10, namely (S s — l)(d n +2 + d n — 2d n +i) = 0. All of the roots A of the 
characteristic polynomial corresponding to this linear relation have modulus |A| = 1, and A = 1 is a triple 
root, which accounts for the growth rate of d n ; the fact that d n — 0{n 2 ) also follows directly from (I3.6p . 
using Y n = 0(1). 

Example 3.7. The tropical version of the recurrence (|2.27l) in Example 12.121 is 

d n+6 + d n = max( 2d n+5 + 2d n+ i , 2d n+i + Ad n+3 + 2d n+2 ). (3.8) 

With the initial conditions d\ = — 1 and d 2 = . . . = cfe = 0, this generates a degree sequence beginning 

-1, 0, 0, 0, 0, 0, 1, 2, 4, 8, 18, 38, 79, 164, 342, 712, 1482, 3084, 6417, 13356, 

that grows exponentially, such that the entropy is £ — log X m ax, where A max ps 2.08 is the root of largest 
magnitude of the polynomial A 4 — A 3 — 2 A 2 — A + 1. To see this, note that the map (|2.28p . in recurrence 
form, is y n +2 Vn = Vn+x + Un+i' so ^ na ^ ^ s tropical analogue is 

Y n+ 2+Y n = \Y n+1 \. (3.9) 
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From the tropical version of Lemma 12.91 an appropriate choice of basis for im B gives the reduction from 
P-8P to (|3.9p , by setting Y n = d n+i — d n+3 — 2d n+2 — d n +i + d n . It can be shown directly that all the orbits 
of (|3.9[) are periodic with period 9 [5], and hence in this case the degrees d n satisfy a linear recurrence of 
order 13, that is (S 9 - l)Y n = (S 9 - 1)(5 4 - S 3 - 2S 2 - S + l)d n = 0. From the periodicity of Y n it is clear 
that d n +A — e?n+3 — 2e?„ + 2 — c? n +i + d n = O(l), which implies that d n ~ CA^ QX for some C > 0. 

Observe that in the two examples above, the tropical recurrences (13. 5p and p.8p both exhibit periodic 
behaviour, in the sense that the maximum on the right hand side is achieved periodically by the first or the 
second entry. If one writes "+" in the case where 5^3=2 [^i>j]+ > S^^I - 6i,j]+ 5 an d "— " otherwise, then 
for the tropical Somos-4 recurrence p.5p one finds that the integer sequence in Example 13.61 repeats the 

block pattern "H 1 1 " of length 8, while for (|3.8|> the repeating symbolic sequence is the block 

"+ H h + H " of length 9. The block length is in accord with the period of the associated periodic 

maps for the variables Y n in each case; for other choices of initial data there are still repeating blocks of the 
appropriate length, but the patterns can be different. Based on a large amount of other numerical evidence, 
we are led to formulate the following. 



Conjecture 3.8. For each recurrence \3. 3]) there exists some k > 1 such that for all real solution sequences 
(d n ), the associated symbolic sequence corresponding to max is eventually periodic with minimal period k. 

Corollary 3.9. For large enough n, every real solution of \3. 3\) satisfies a linear recurrence relation of order 
kN with constant coefficients. 



The above corollary results from the fact that, given Conjecture 13. 8[ for large enough n one can regard 
(13.3[) as being equivalent to the iteration of a linear recurrence relation whose coefficients vary with period 
k, and then the following result can be applied. 

Lemma 3.10. Suppose that a sequence (s„) satisfies a linear recurrence of order i whose coefficients are 
periodic of period k, say 

e-i 

s ±, = V c (r) s ^ c (r) - c (r) 

r=0 

Then the terms of the sequence also satisfy a linear recurrence of order ki with constant coefficients, of the 
form 

e-i 

Sn+M — ^ Cr S n + r k- (3.10) 
r=0 

Proof: This result should be well known in the literature on linear recurrences, but for completeness we 
sketch the proof. It suffices to consider the (£ + 1) x (£ + 1) matrix 



( S n S n+ i . . . S n+ l \ 

Sn+k Sn+fc+1 ■ • • $7i+k+£ 

\ s n+ ke Skt+i ■ ■ ■ s„+(fc+i)^ / 



(3-11) 



which has vanishing determinant, by virtue of the fact that the vector (c4 , • ■ • , c« l \ — 1) T is in the right 
kernel. The recurrence with constant coefficients is obtained from the left kernel (the kernel of $^). □ 



In general, it is not easy to determine the period of the symbolic sequence, and Conjecture 13.81 seems 
challenging. The real orbits of the recurrences p.3p can be quite complicated. 
A further refinement of Conjecture 13.81 is possible. Let the matrices 
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be defined in terms of the palindromic vectors a± of size N — 1, which correspond to the exponents [±Oj]+ = 
[±&ij_l_i]+ (where T is a zero row vector of size iV — 1, and ljy— 1 is the corresponding identity matrix), so 
that their characteristic polynomials are P± as in (|3.1[) . respectively. Also, let 

n= f[ M M. ...M, 

j=i,...,k 

be the path-ordered product corresponding to the symbolic sequence of length k for a particular orbit of 
(|3.3[) , defined by an appropriate choice of e.j — ±, and let p(M±) and p{U) denote the corresponding spectral 
radii. 

Conjecture 3.11. Let k be the period of Conjecture \3.8\ There are two possibilities: 

(1) // |A maa; | = p(M.±) > /9(Mzp) > 1 then k = I, with repeated symbol e\ — ±, respectively. 

(2) 7/|A max | = p(M+) = p(M_) then p{U) = |A max | fc with k > 1. 

To see how Conjecture 13.11 now follows from all the rest, consider the matrix in the proof of Lemma 
13.101 for the case s n = d n and £ = N. If $„ is the N x N submatrix such that det <fr n is the upper left 
connected minor of size N in (|3.11l) . then a single iteration of (|3.3p gives = <&„]V[ ei for some choice of 

e% = ±. After k iterations it follows from Conjecture 13.81 that $„+fe = $ n II, and if n is large enough then 
'i'n+rfe = $nll r for all 7*, for some block of symbols £j of length fc that is fixed up to a cyclic permutation 
(which depends on the choice of n mod k). By the Cayley-Hamilton theorem, if P(k) = k n — ^2^Sq c r n r is the 
characteristic polynomial of II (which is independent of cyclic permutations of the M e . ), then P( l S fc )$„ = 
<f>„P(II) = 0, which shows that $„, and hence d n , satisfies the recurrence (|3.10[) with I — N, for all n 
large enough. Thus the characteristic roots A for the growth of d n satisfy P(X k ) = 0, which implies that 
d n ~ C|A moa; | n for |A moa: | = p(II) 1 / fc > 1, or d n has polynomial growth for |A mQ2; | = 1. If k = 1 holds, then 
in either case (1) or case (2) of Conjecture 13.111 II = M + or M_, so |A maa; | is given by one of the spectral 
radii of the matrices M±, whichever is the larger, and £ = log |A maa; |. The condition in case (2) is required 
to ensure that first part of the statement of Conjecture 13.11 holds for k > 1 as well. 

For the second statement in Conjecture 13. 1[ note that P± in (|3.1j) are reciprocal polynomials (their 
coefficients are palindromic), so that P+(A _1 ) = \~ N P+(A), and similarly for P_. Let 

N 

s± = Xl[ ±5l J+ 1 ]+' 
J'=l 

respectively. By the symmetry of the matrices B — > —B (or equivalently, the freedom to replace the quiver 
Q by its opposite), it can be assumed without loss of generality that S+ is the larger of the two, so S+ > SL, 
and take S+ > 3 so that condition (|3.2I) holds. Now P + (0) = 1, and P+ (1) = 2 — 5+ < —1, so P + has a real 
root between and 1. The reciprocal property implies that the reciprocal of a root is also a root, hence P+ 
has a real root larger than 1, implying that |A maa; | > 1 and the entropy is positive. 

The cases for which (I3.2[) does not hold are easily enumerated, and it can be checked directly that 
|A ma o;| = 1 m these cases. The main conclusion of this section is the following. 

Theorem 3.12. Suppose that Coniecture \3.1\ holds. Then up to the symmetry «S+ -h- S-, there are only four 
distinct choices of the pair (S+,S-) for which the algebraic entropy is zero, corresponding to four different 
families of recurrences: 

(i) (S+, S-) = (1, 0): For even N — 2m only, the recurrence is 

*£n+2m -^n ^n+m ~\~ 1* (^'12) 

(ii) (S+,S-) = (2,0): For each N > 2 and 1 < q < [N/2\, the recurrence is 

X n +N X n = X n +N-q %n+q + 1- (3.13) 
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(iii) (£+, S-) = (2, 1); For even N = 2m onfa/, and 1 < g < m — 1 ? i/ie recurrence is 

•£n+2m ~~ *£n+2m — g ~l~ ^n+^' (3.14) 

(iv) (S+, = (2, 2); For each N > 2 and 1 < p < q < [N/2\ ? the recurrence is 

(3.15) 

The simplest is case (i), where the recurrence (|3. 121) decouples into m independent copies of the Lyness 
map: all the orbits are periodic, and the overall period of the sequence of x n is 5m. For each n the function 
F n = x n + x n+m + x n+ 2m + ^n+3m + in+4m is invariant, and can be written as a function of x„ and 
x n+m only, using (I3.12[) . Moreover, the functions F\, F2, ■ ■ ■ , F m are independent and Poisson commute with 
respect to the bracket corresponding to the symplectic form w, so trivially this system is also integrable in 
the Liouville- Arnold sense. 

The integrability of the families (ii),(iii) and (iv) above is discussed in subsequent sections. 



4 Linearisable recurrences from primitives 

The primitives introduced in |13j are the simplest examples of cluster mutation-periodic quivers with period 
1, and are the building blocks of all such quivers. The nonlinear recurrences that arise from the primitives 
have the form (|3.13l) , corresponding to case (ii) in Theorem 13.121 and can be rewritten as 

Xn+N %n = %n+p %n+q + 1, P + q = N, (4.1) 

so that for each q = 1, . . . , [N/2\ there is a different recurrence corresponding to the primitive Pff . When 
p and q are coprime, the associated quivers are mutation-equivalent to the affine Dynkin quivers A q ^ p , 
corresponding to the ^.^Li Dynkin diagram with q edges oriented clockwise and p oriented anticlockwise, 
while if gcd(p, q) — k > 1 so that p = hp, q = kq, then the quiver is just the disjoint union of k copies of 
Ag } p. In the latter case, it is clear that the recurrence (|4.ip is also equivalent to k copies of the recurrence of 
order N/k corresponding to the coprime integers p, q. Hence it is sufficient to consider only the case where 
p, q are coprime, which we will assume from now on. 

The cluster algebras generated by affine A-type Dynkin quivers arise from surfaces, and are of finite 
mutation type [5], meaning that only a finite number of distinct quivers is obtained under sequences of 
mutations (|2.1j) . However, by the classification result in [TU], these cluster algebras are not themselves 
finite: the recurrence (|4.1I) generates infinitely many cluster variables starting from the initial cluster x = 
(xi, . . . ,xjv). 

It was conjectured recently that the cluster variables in cluster algebras obtained from affine Dynkin 
quivers satisfy linear recurrence relations, and this was proved for all but the exceptional types pQ. A 
different proof using cluster categories, valid for all affine Dynkin types, was subsequently found by Keller 
and Scherotzke [53] . In the case of A q<p (with coprime p, q) a proof of the corresponding linear recurrence 
relations was already given in |13) . Here we present a more direct derivation of the linear recurrences arising 
from these affine A q p quivers, before using the Poisson structures from section 2 to explain the Liouville 
integrability of the maps defined by (14.11) . 

4.1 Linear relations with periodic coefficients 

The key to the properties of the nonlinear recurrence (I4.1[) is the fact it can be written in the form 

det 9„ = 1, where = ( Xn \ . (4.2) 

\ %n+p X n+N J 

The identity (|4.2|) is the frieze relation (see e.g. [T]); it implies that the iterates of (|4.1j) form an infinite 
frieze. 
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Upon forming the matrix 



Xn X n -i-q X n +2q 

&n = | x n+p %n+N Xn+N+q 
X n +2p ^n+JV+p %n+2N 



(4.3) 



one can use the Dodgson condensation method [6] to expand the 3x3 determinant in terms of its 2 x 2 
minors, as 

det §„ = — !— (det det ® n+N - det * n+g det *„ +p ) = 0, 

X n +N 

by (|4.2[) . By considering the right and left kernels of we are led to the following result. 
Lemma 4.1. The iterates of the recurrence \4-l\ ) satisfy the linear relations 

X n +2q - Jn X n +q + X n — 0, (4.4) 
X n +2p - K n X n+p + X n = 0, (4.5) 

whose coefficients are periodic functions of period p, q respectively, that is 

J n +p = Jn, K n+q = K n , for all n. 

Proof: A vector in the (right) kernel of can be written as k„ = (J n , — J n , 1) T , where the third entry is 
normalised to 1 without loss of generality. From the first two rows of the equation ^„k„ = we have the 
linear system 

— Jn \ _ ( X n+ 2q 
Jn J V X n+N+q 



From Cramer's rule it follows that 

J n = -(det*„)" 1 

and 



Xn+2q X n -\-q 
Xn+N+q %n+N 



(det# n ) x det* n+9 = 1, 



J n 



Xn X n J r 2q 
Xn+p *£n+7V+q 



where we made use of (|4.2I) . Hence the recurrence (|4.4p is given by the first row of ^>nkn = 0. The second 
and third rows of the latter equation provide the 2x2 linear system 



+p 



J n 
J n 



Xn+N+q 
X n +2N 



whose solution gives an alternative formula for J„. Noting that the coefficients in this linear system are the 
same as those in the first two rows of the equation ^ n+p k n+ p = 0, it follows that J„+p = J n for all n. If we 
consider the left kernel of \P„ (the kernel of ), then by symmetry we obtain the recurrence (|4.5[) . whose 
coefficient K n has period q. □ 

Remark 4.2. In the particular case q — 1, corresponding to Ai t N-i Dynkin quivers, the coefficient K n has 
period 1, so K n+ \ = K n = JC (constant) for all n. In that case the recurrence (|4.5[) becomes 



X n +2N-2 + X n — JC X^+jv-I) 



(4.6) 



which is a constant coefficient, linear difference equation for x n . This was first shown in [13] and is precisely 
of the form derived for affine A-type quivers in [TJ [23] . 
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The quantity K. in (14.61) can be considered as a Laurent polynomial of the initial data Xi, X2, ■ ■ ■ , xjy for 
the map ip associated with (|4.1| ; from (|4.6p this is obtained explicitly by repeatedly using the nonlinear 
recurrence in order to write K in terms of lower order iterates until it is given in terms of the N initial data. 
Moreover, the fact that it is independent of n means that K, is a first integral of ip. 

In the case of general coprime p, q, it is straightforward to obtain first integrals of the map <p corresponding 
to (|4.1|) . Indeed, the orbit of J\ under the action of ip generates p functions J±, J2, ■ ■ ■ , Jp, which can be 
written as rational functions of the N initial data via back-substitution in the relation 

J n = X " +X " +2g . (4.7) 

Clearly any cyclically symmetric function of J\, J2, . . . , J p is a first integral of tp and the first p elementary 
symmetric functions of these variables are independent functions of the J n . Given that Ji, J2, . . . , J p are 
functionally independent as functions oi x±,X2, ■ ■ ■ ,xn, any p independent cyclically symmetric functions of 
the J n can be picked as first integrals. Similarly, the action of tp generates the functions K\, K2, . . . , K q , and 
any q independent cyclically symmetric functions of them provide another set of independent first integrals of 
(p. This gives a total of p + q = N first integrals. These cannot all be functionally independent, because that 
would imply that the map is periodic. However, the generic orbit generated by (|4.ip is not periodic. Below 
we describe a single functional relation between these first integrals (equation (|4.12j) below), by considering 
the monodromy properties of the linear relations (|4.4[) and (|4.5|) . 



Remark 4.3. The existence of N — 1 independent first integrals Ij, together with the volume form Q in 
(|2.17p . means that the map ip is Liouville integrable in a rather elementary sense. By applying a method in 
[4], one can pick, say, the first N — 2 integrals and obtain a Poisson bivector field J which is invariant (or 
anti- invariant, for odd N), namely 

J = rijdlij idlAT_ 2 , 

where the N- vector field il is defined by = 1. By construction, this Poisson structure has Casimirs Ij 
for j — 1, . . . , N — 2, and restricting to the two-dimensional symplectic leaves one has an integrable system 
with one degree of freedom, with Ijv-i being the independent first integral. However, we shall see that the 
Poisson structure coming from the two-form (|2.10|) leads to more interesting integrable systems. 

4.2 Monodromy matrices and linear relations with constant coefficients 

The relation (|4.4j) implies that the matrix \l/„ satisfies 

' -1 



y n+q =$„L„, L„ j y (4.8) 

Upon taking the ordered product of the L„ over p steps, shifting by q each time, we have the monodromy 
matrix 

M„ := L„L„ +( j . . . L n+ ( p _ 1 ) 9 = ^n+pq- (4.9) 
On the other hand, the recurrence (|4.5j) yields 



<f„ +p = L„ L n = ( _° x ^ J , (4.10) 

which gives another monodromy matrix 

M„ := L n+ ( g _i) p . . . L„ +p L„ = & n +pq ■ (4-11) 

The cyclic property of the trace implies that 

K n := trM„ = trM n . (4.12) 
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Also, since L n has period p, shifting n — > ri + p in (|4.9p and taking the trace implies that JC n + p = }C n - 
Similarly, from (|4.11[) we have lC n + q = fC n . Now the periods p and q are coprime, and since fC n has both 
these periods it follows that /C„ = K, ^constant, for all n, hence K. is an invariant of ip. 

From the expression (|4.9[) it further follows that K, is a cyclically symmetric function of the J n , n = 
1, . . . ,p, while from (|4. 1 1[) it is also a cyclically symmetric function of the K n , n = 1, . . . , q. Thus we see 
that the equality of the traces in (|4.12l) provides the aforementioned functional relation between these two 
sets of functions. 

We are now ready to show that in general the iterates of (I4.1j) satisfy a linear relation with constant 
coefficients. The existence of such a relation follows immediately from (|4.4[) or (|4.5p . by applying Lemma 
13.101 but the monodromy matrices provide more detailed information about the coefficients. 

Theorem 4.4. The iterates of the nonlinear recurrence j4.1\ l satisfy the linear relation 

Xn J \-2pq X n ^ ^n+pep (4T3) 

where K, is the first integral defined by with K. n = K. for all n. 

Proof: Using (|4.9I) we see that v f n +p q — ^nM„, so \E I ' n +2p q — ^ n^-n^-n+pq — ^ n^-n by periodicity. Noting 
that M„ is a 2 x 2 matrix, with det M„ = 1 and trM n = /C, yields (by Cayley-Hamilton) 

^n+2 Pq - K.^n+p q + *n = $n(M 2 n - JCM n + I) = 0. 

The (1,1) component of this equation is just (|4.13[) . □ 

By making use of the Chebyshev polynomials of the first and second kind, defined by T n (Q = cos(n9) and 
U n (Q = sin((n + 1)9)/ sin6>, respectively, for C = cos0, the linear equation (|4.13p yields an exact expression 
for the iterates of the nonlinear recurrence. 

Corollary 4.5 (Chebyshev polynomials). The recurrence has the explicit solution 

Xi+n P q = XiT n (K,/2) + ( x i+pq 1 —\ U n -i(K,/2), i=l,...,pq-l, for all n. (4.14) 

4.3 The structure of monodromy matrices 

Here we give some relations between the elements of M„ , which provide properties of a natural Poisson tensor 
associated with the functions Ji. Analogous results regarding M„ and Ki also hold, since the structure of M„ 
and M„ is the same up to switching p «-» q, Ji «-» Ki and taking the transpose. To simplify the presentation, 
we concentrate on the case of . (Similar results hold for P 2 mi witli p + q = 2m and p, q coprime.) The 
remarkable fact is that (when is even) these properties of the Poisson bracket are derived directly from 
the monodromy matrix. In subsection 14.41 we proceed to show how the Poisson algebra of the Ji is derived 
from the Poisson bracket between the coordinates Xj . 

For the case q — 1, we denote the matrix M„ = L„L„ +1 • • • L n+p _ 1 by 'M. n 2m \ when p = 2m — 1, m > 1, 
and Mn + , when p = 2m, m > 1. 

It is important in the calculations below that M r ^f +1 ^ depends only upon the variables J„, . . . , J n+p _i. 
For the moment, this is not really a monodromy matrix, if we do not assume any periodicity. The calculations 
below give a recursive procedure for building the matrices M^ 2m ^ and M^ 2m+1 ^. 

The recursion Ml P+3) = m1 p+1) L ii+p L„ +p+ i, with the short-hand notation A = A^ p+1) , A = Al P+3) , 
etc., leads to 

A = J„+ p B — A, B = (Jn+pJn+p+l — 1)B — J n +p+iA, 

(4-15) 

C = — C, D — (Jn+pJn+p+l — 1)D — J„+ p +iC, 

so that K- = —K. + J„+pB — Jn+p+iC + J n + p J n + p +iD. 
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Lemma 4.6 (Relations for M). The components o/m! p+1 ' 1 satisfy the relations 

a An = gg„ = Q A (p+1) _ B (p+1) = ^ C (p+D _ C ^ 

where fC^ = Al P+1) + Dl P+1) . 

(2) 

Proof: We just prove this for the even case. For m = 1, M„ = L„, which clearly satisfies these relations. 
The recursion f|4.15[) provides us with an inductive step. We have 

dk 8B dA dB dB 

Tr-T~ — Jn+2m~lTT-; Try — U, 7—— — (Jn+2m- Hn+2m — 1)^77 — I" Jn+amTTT" — U. 

Then 

<9C 9/C dC ~ 

Tpr" — Jn+2m-lTr? Try- — A ~ Jn+2m-\B — — A, 



dU _ d(A+D) _ dK_ QT1 ^ d£ 

Finally 



where we have used jj£ = a ^J^> = jjjfc, and dJ ^ n = -C + J n+2m _iD = C. 



dfC _ _dK_ j <!£_ ^Sl 

8J n ~ ~ 0J n + Jn+2m - X dJ n ~ Jn+2m dJ n + J n+2m-lJn+2m ^ 

= (•4+2m-l^ij+2m — l)"^ - ? ^ ^n+2mA = — ( •/i ! +2m-l'4+2m — 1)B + J„-|_2mA = — B, 

OJn 

where again we have used -gj- = -^j- . □ 

We can use the above relations to form recursion operators which can be used to build the functions 
K.( p+1 \ Starting with K.^ = J n and K,^ = J n J n +i - 2, we can use /C(f+3) = Tl^K^ +1 \ where the 
recursion operator is 

d 2 d d 

TZ(P) = J n+p J n+p+ i Jn+p-yj Jn+p+lJTJ h (Jn+pJn+p+l — 1)- (4.16) 

OJ n OJ n +p-i OJ n OJn+p-1 

Remark 4.7. An alternative formula for K, is given by a link with the dressing chain: 

P / f)2 \ P 

^n(-^)rj>- 

When p is odd, this formula follows from the results in [35], by setting — >■ and gi — > Ji. 



4.3.1 Link with the Poisson structure 

In the case TV = 2m, we can use the monodromy matrix to build the Poisson bracket for the functions J„, 
Staying within the context of one quiver P 2 to f° r fixed m, we now reinstate the periodicity L n+p = L„. 
From (|4.9I) . by periodicity, for all k and n we have 

M„ + feL„ + fe = L n _|_fcM n _|_/._|_i. 

When k — 0, this implies 

Ajj+i = D„ — J„C„+i, B„ + i + C n = J„(D„ — D Jl+ i), C„+i = — B„, D„ + i = A„ — J„B„. 
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Shifting n, we can write C n = — B„_i = — B„ +p _i. Only two of the remaining equations are independent, 
leading to 

JrJin = A n — D„ + i, B„ +p _i — B„ + i = J„(D n+ i — D„). 

The equations for k ^ are obtained by shifting the indices. The first equation leads to 

p— 1 p— 1 

(~ l) fc+1 J n -|_fcB n+ fc = l) fc+1 (A„ + fc — D„ + fc + i) = A n+ x — A n+P = A„ + i — A„, 
fe=i fe=i 

since the remaining sum consists of 

m— 1 m — 1 

(A„+2 3 +i + D n+2 j+i) - (A„ +2 j + D«+2i) = ("i - 1)£ - (to - 1)/C = 0. 
i=i 3=1 

On the other hand B„ +p _i — B„ + i = J„(D rl+ i — D„), so 

p-i 

>/n ^^{ — ^-) k+1 Jn+k^n+k + B Jl+p _i — B Jl+ i = J„(A„ + i — A n+p + D„ + i — D n ) = 0. 
k=l 

If we take cyclic permutations, we obtain a matrix equation of the form 

(p(2) +P (o)) b = 0) with b=(B„,B„ +1 ,...,B„ +p _ 1 ) T = -V/C, 

where the components of the matrix P^ 2 ' + Pw are those of the Poisson tensor in Lemma T4.9I below, and 
we have used the formula B n = —■§j L - for all n, coming from Lemma 14.61 
We can summarise these results in 

Theorem 4.8. The function JC, defined by fi4-12\ ) is the Casimir of the Poisson bracket of Lemma \4 ■ 9\ 

(p( 2 )+p(°))V/C = 0. 

4.4 Poisson brackets and Liouville integrability for P 2 ^ 

It was proved in |14] that the linearisable maps coming from the primitives P^ (the A\.n-i Dynkin quivers) 
are Liouville integrable when N is even. We give the proof here, since it is the basis for understanding the 
integrability of the maps for the other Pjp quivers. 

The proof starts from the Poisson bracket for the cluster variables. The matrix B in this case is nondc- 
generate, having the form 

N-l 

B = t n - tJj, with t n E r+l,r - Ei.jv, (4.17) 

r=l 

where E r:S denotes an element of the standard basis for gl(N). The "skew rotation" matrix rjv plays an 
important role in the classification presented in [15] . The matrix (|4.17j) for N — 4 is obtained by setting 
c = in igUH) . 

By Theorem 12. 6i the map is symplectic, and hence there is a nondegenerate Poisson bracket of the form 
(|2.6|) . with C = B^ 1 , up to scaling. In accordance with [14], we take 

f N-2s+l 

C = Tjy + (t^) 3 + • • • + {tn) N 1 = ^] J]] (E rjr +2s-l — E r +2s-l,r) j 

s— 1 r— 1 
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so that CB = 21, which gives 

f sgn(fc- j)xjX k , k-j odd, 

{x j ,x k } = \ (4.18) 
I otherwise, 

for j, k = 1, . . . , N, with sgn denoting the sign function. 

The key to the Liouville integrability of the P^ maps is the expression of the Poisson bracket between 
the periodic functions J„, which appeared in [21] (for N = 4) and [14] (for general even N = 2m). 

Lemma 4.9. For even N = p + 1 and q = 1, the functions J n given by |f^.7p define a Poisson subalgebra of 
codimension one in the algebra with brackets 

{ Jj, J fe } = 2sgn(fc - j) (-iy+ k + 1 Jj J k + 2(6 jik+1 - 8 j+1>k + 8 j+p -im - fy.k+p-i). (4.19) 

for j,k = l,...,p. 

The bracket for the J n is clearly a sum {,} = {, } 2 + { , }o, corresponding to the splitting of the Poisson 
tensor into homogeneous parts P^ 2 ) +p(°), as in Theorem 14.81 The quadratic bracket { , }2 is log-canonical, 
while the degree zero bracket can be defined simply by specifying its only non-zero terms as 

{ J j+1 , J, } = 2 = -{ J h J j+1 } (4.20) 

for all j mod p = 2m — 1. Since (14. 19)) is still a Poisson bracket after scaling Jj —> /iJj for arbitrary /i, the 
brackets {, }o and {, }2 are compatible, so define a bi-Hamiltonian structure. This means that one can use 
the standard bi-Hamiltonian chain [26 , defining a sequence of functions Ij which satisfy 

{Zj,?k}o = = {lj,l k } 2 , for all j, k, 

where the sequence starts from Iq = • Jj, the Casimir of the bracket {, }o given by (|4.20l) . and finishes 
with Im-x = Y\j Jj, the Casimir of {, By Theorem l4.8[ the function JC defined by (|4.12l) is the generating 
function for these integrals, so that 

m— 1 
3=0 

where Tj is the term of degree 2j + 1 in the variables Ji . Since these integrals commute with respect to both 
brackets, they commute with respect to the sum { , } 2 + { , }o, and hence provide m commuting integrals 
for the map ip of the variables Xj in even dimension N, which implies that the map is Liouville integrable. 
In summary, we have 

Theorem 4.10. For N = 2m and q = 1 the map ip defined by has m functionally independent Poisson 
commuting integrals, given by the terms of each odd homogeneous degree in the quantity K,, as given by 
equation {^.21^ . The map is also superintegrable, having a total of N — 1 independent first integrals. 

As discussed in subsection 14.21 extra first integrals are obtained by choosing m — 1 additional cyclically 
symmetric functions of J\ , . . . , J p . 

4.5 Primitives of the form P 2 m+i 

The recurrences (|4.1j) for p = 2m and q=l are given by 

X n +2m+l X n = X n+ 2m X n +1 + L (4.22) 

The formula (|4.17l) still holds, but now the matrix B is singular. It has a one-dimensional kernel, spanned 
by the vector u = (1, — 1,1) T , which generates the scaling symmetry 



X ( - iy+1 Xj, AeC*, (4.23) 



21 



and im B is spanned by 

v j =e j +e j+l , j = l,...,2m, (4.24) 
where is the jth standard basis vector. Hence, by Lemma \2. 71 the coordinates 

yj = XjX j+ i, j = 1, ...,2m, (4.25) 

are invariant under the scaling (|4.23p . and the degenerate form (|2.10l) pushes forward to a symplectic form 
(I2.22[) in dimension 2m, whose coefficients bjk are the matrix elements of 

p> _ 2 | 3 _ , „2m-l 

^ '2m '2m "T T 2m ' ' ' ' T 2m > 

where T2 m is the 2m x 2m version of t^. The inverse of this is the skew-symmetric matrix 

B- 1 = rj m + (rj m ) 2 + (rj m f + • • • + (r^) 2 - 1 , 

with all components above the diagonal equal to 1, giving a nondegenerate Poisson bracket for the yj, i.e. 

{yj,Vk} = sgn(fc- j)yjVk, l<j, fc<2m. (4.26) 

Upon applying the rest of Theorem 12.61 we see that (|4.22l) induces a symplectic map on the variables yt , 
given (for m > 2) by 

■ (2/l,y2,.-.,y2m-l,2/2m) l"> (?/2, 2/3, • • • , 2/2m, */2ro+l), ( 4 - 27 ) 

where 

2/2m+l = 2/2J/4 ' ' ' y2m(y22/4 ' ' ' 2/2m + 2/3 ' ' ' V2m-l) / {VlVl ' ' ' 2/lm-l)- 

The map is simpler for the case m = 1, given in Example 14.111 below . 

By the general discussion above, the iterates x n satisfy the linear relation (|4.6[) . where K is the trace of 
the monodromy matrix M„. The latter is given in terms of the quantities 

J n = Xn+Xn+ \ n = l 2m, (4.28) 

which cycle with period 2m under the action of the recurrence (14.221) . The polynomial K, can be expanded 
as 

m 

fC = ^2(-l) m+j l j , where X = 2, (4.29) 

3=0 

and each polynomial Ij is homogeneous of degree 2j in the variables J n . The non-trivial homogeneous 
components I\, . . . ,I m provide m first integrals for (I4.22[) . and an additional m independent first integrals 
can be obtained by choosing cyclically symmetric functions of J„ for each odd degree 1,3,..., 2m — 1. 

However, not all of these first integrals reduce to functions on the 2m-dimensional symplectic manifold 
with coordinates yj . Note that the scaling symmetry (|4.23[) acts on the variables (14.28[) according to 

Jn — > X 2 ^ n+1 J n . (4.30) 

Applying this to the formula (|4.16l) shows that each component of K, is invariant under scaling, which means 
that the m first integrals X\ , . . . , T m can be rewritten as functions of the scale- invariant variables yj . The 
Liouville integrability of the map (I4.27P then follows, provided that these m functions are in involution with 
respect to the bracket (|4.26[) . We shall not pursue the case of general m further here, but content ourselves 
with presenting some low-dimensional examples. 
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Example 4.11 (The primitive P 3 )• For m — 1 the map tp for the Xj variables is (|2.8p , which is associated 
with the affine Dynkin quiver. The matrix B has rank two, so in terms of the variables yj given 
by (|4.25l) the symplectic form has the log-canonical form (I2.26|) . and the induced map of the plane is 

tp ■ ( Vl ^ h> f 2/2 

Symmetric functions of the period 2 quantities 

_ Xl + X 3 XiX 2 + x 2 x 3 + 1 

J\ — , -h = 

X 2 XiX 3 

give two first integrals for the map (|2.8p . namely 

J 1 + J 2 and X 1 = J 1 J, = K + 2= {vi+y2){yi+y2 + 1) . 

2/12/2 

The latter is defined on the (2/1,2/2) plane, and <p*I\ = X\, so the symplectic map tp with one degree of 
freedom has an invariant function and hence is integrable. 

Example 4.12 (The primitive -P5 )• For m — 2 the recurrence (|4.22[) has the first integrals Ii and I2 given 
by the homogeneous terms of degree 2 and 4, respectively, in the expression (|4.29[) : 

K = 2 - ( Ji J 2 + J 2 J 3 + J 3 J 4 + J4J1) + J1J2J3J4 = X -X 1 +X i , (4.31) 

where the J„ are defined by (|4.28j) . Picking another pair of cyclically symmetric functions of J n , of degrees 
1 and 3, say, adds two more independent first integrals. 

Now defining yj for j = 1, 2, 3, 4 by (|4.25p . these variables are endowed with the nondegenerate Poisson 
bracket (|4.26p . which is invariant under the map (j4,27[) : 

0- (2/1,2/2,2/3,2/4)^ (2/2, 2/3, 2/4, 2/22/4(2/22/4 + 2/3)/(2/l2/3))- 

To show that this map is Liouville integrable, it is necessary to verify that {Zi,^} = with respect to this 
bracket. 

The terms in the formula (|4.31l) can all be expressed via the functions 

Wi = JiJ i+1 , (4.32) 

so that Zi = wi + W2 + W3 + 11)4, I 2 = W1W3. From (|4.30p it is clear that these Wi are invariant under the 
action of the scaling symmetry (|4.23p . This means that they can be written as functions of yj, viz: 

(2/1 +2/2) (2/2 + 2/3) (2/2 + 2/3)(2/3+2/4) (2/3 + 2/4) (2/2 2/3 + 2/1 2/3 +2/ 2 2/4) 

Wl = 5 , W 2 = 2 > W 3 = 2 • 

2/2 2/3 2/12/32/4 

Under the action of tp, since the J n cycle with period 4 under tp, the transform as 

-» .» W1W3 

tp Wi = U>2, tp W2 — W3, tp W3 = W4 = , W4 — W\. 

Although only the first three are independent, it is convenient to make use of W4 as well. 

The first three Wi form a three-dimensional Poisson subalgebra of the yj, which is non-polynomial: 

r l r l WlW3 r 1 

\WX,W2} = W1W2 - Wl - W 2 , \Wi,Wz\ = W2 , {w 2 , W 3 } = W 2 W 3 - W 2 ~ W 3 . 

W2 

The Casimir of this algebra is Xi — X% = Wi + w 2 + w 3 + w ^ 3 — wiw 3 = 2 — K. 

Since Xi and X 2 are both functions defined on this subalgebra, it follows that {Xi,X2} = {Ii,K,} = 0, so 
the two first integrals are in involution, as required. 
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Remark 4.13. The Poisson bracket of the four functions u>i can be calculated in polynomial form as 

{Wl, W 2 } = WiW 2 - Wi - W 2 , {wi,W 3 } — W 2 — W4, {wi,U>4,} = tOi +W4, — WiWi, 

with the remaining brackets following from the cyclic property. This bracket has the two Casimirs 

C\ = W1W3 — w 2 W4 and C 2 = w\ + w 2 + W3 + W4 — W1W3, 
so that the three-dimensional algebra for w±,w 2 , W3 arises from the constraint C\ = 0. 

4.6 Primitives Pjfi with q > 1 

For general q, to obtain the quiver Pfi' we must modify the formula (|4.17|) and write 

B = r q N {rl) T . (4.33) 

In this case, we have to take into account both sets of functions J„, n = 1, . . . ,p and K n , n = 1, . . . , q, 
defined through (|4.4[) and (|4.5p respectively, with the function K. being given by the two formulae of (|4.12[) . 
It turns out that the essential properties of the quantities J n with period p can be obtained by considering 
the J n in the case of Pp+i (the Ax tP quiver), and applying a suitable permutation of indices; and similarly 

the properties of the K n are the same as those of the J n for Pg+i, up to a permutation of indices. 

Concentrating for the moment on the functions J n , we consider the formula (|4.9j) for M n . Each of the 
matrices ~L n +i q (after using the cyclic property J n +p = Jn) is j us t one of the matrices L„, n = 1, . . . ,p. For 
coprime p, q, each of the J n appears exactly once in this product and in a specific order, which defines a 
permutation a of the integers (1, . . . ,p). For this discussion, for a general pair p, q (with q < p) let us use 
lC Pyq to mean tr M n , considered as a function of the J n . Then we have 

= fcp,q(Jl,- ••,Jp) = £p,l(J<7(l)> • • • , J<t( p ))- (4.34) 

Similarly, in terms of the functions K n , the formula (I4.11[) for M„ defines a permutation a of (1, . . . , q), and 
if we write /C 9lP to denote trM„, then we have 

K = )Cq tP (Ki, . . . , K q ) = /C g> i (#£(!),..., K&( q )). (4.35) 

(There is no risk of confusion between /C 9 , p and K- Py q once we have fixed q < p.) 

For more detailed properties of the maps defined by (|4.1[) . and the associated quantities J n and K n , it is 
necessary to consider even/odd N separately. 



4.6.1 The even case 

In the case that N is even, when q and p are coprime it can be shown by elementary row/column operations 
on the matrix (|4.33|) that det B — 4. The matrix B is of Toeplitz type, and invertible, with an inverse of 
the same type that defines a nondegenerate Poisson structure of the form (|2.6|) . With the choice of scale 
C = 2£?~ 1 , the Poisson bracket for the Xj is given explicitly by 

( Bga(k-j)(-l) r+s XjX k , k-j odd, 

{Xj,x k } = < (4.36) 
I otherwise, 

where (by the Euclidean algorithm) the integers r, s are uniquely determined by writing 

— (p — \k — j\j = sm — r£, for < r < m, < s < £, 
in terms of the coprime integers £ — (p — q)/2, m = (p + q)/2. (Note that p and q are both odd.) 
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Consider the functions J n once again, and for fixed n define the following sequences: 

Y — ■ — 10 7t_T _ + X k+2 , _ 

Aj — x n+(j—l)qi J — l j A • ■ • i J k ~ J n+(k-l)q — > K — L i *"> ■ ■ ■ ■ 

Note that the sequence of j\ is also periodic with the same period: j\ +p = J k , in fact (up to the choice of n, 
which just gives an overall cyclic permutation) the ordering of this sequence corresponds to the permutation 
a in (|4.34j) . i.e. j\ = J a / k )- Using the Poisson bracket (|4.36j) . we compute {X\,X2~\ = {x n ,x n +q\ = XyX^, 
{Xi, X$} — {x n , x n +2q\ — 0, and so on, and then in each case we find that, for a suitable range of indices, 
the Xj satisfy the same Poisson algebra (|4. 18[) as the Xj in the 

case, so that the bracket is 



{ Xj , X k } 



"3 " x ± p+i 
sgn(fc- j)XjX k , k~j odd, 

otherwise, for \j — k\ < p. 



To verify the analogue of Lemma l4~9) it is sufficient to calculate the brackets { j\, j|} for j — 2, . . . , (p+ 1 )/2 
and then use the periodicity, which shows that the permuted quantities j\, satisfy 

{ J], 4 } = 2 sgn(fc - j) {-\y +k+1 j]jl + 2(6j, k+1 - S j+llk + 5 j+p ^ llk - ^, k+p-i), 

for j, k = 1, ... ,p. This is the same as the Poisson algebra (|4.19j) for the J„ in the case of -Pp+i< An identical 
argument implies that, up to the permutation a, the q functions K n satisfy the algebra (|4.19p corresponding 
t0 P (1) 

Remark 4.14. Theorem 14. 8[ together with (|4.12|) . implies that the same function JC is simultaneously the 
Casimir of the J n subalgebra and of the K n subalgebra. 

By further direct calculation using the bracket (|4.36p one can verify that 

{Ji,Ki} = Q Mi =^ {Ji,Kj} = Vi,j, (4.37) 

where the second statement follows from the first by repeatedly shifting i i + 1, and using periodicity and 
the coprimality of p and q. 

Thus, from (I4.37[) , we see that these two subalgebras Poisson commute with each other. The J„ subalgebra 
provides (p+l)/2 commuting integrals, which are found by applying the permutation a to each homogeneous 
component of the sum (|4.2ip for the case of the primitive P^+i ■ Similarly, by applying the permutation a one 
obtains (q + 1)/2 commuting integrals for the K n subalgebra. These two sets of integrals also commute with 
each other, and the relation (|4.12[) provides a single constraint, which gives a total of (p+l)/2+(g+l)/2 — 1 = 
N/2 independent commuting integrals, as required for Liouville integrability. 

Example 4.15 (The primitive P^)- The matrix C = 2£>~ 1 is Toeplitz, with top row 

(c u ) = (0,-1, 0,1, 0,1, 0,-1), 

which specifies the Poisson bracket (|4.36|) in this case. This Poisson bracket is invariant under the map 

tp: (xi,...,xs) i-> (x 2 ,...,£ 9 ), xg = XiX6 + 1 . (4.38) 

x-y 

The functions J„, which cycle with period 5, are given by 

Xy + Xj X 2 + Xs T XyX 3 + X4XQ + 1 X2X4 + X 5 Xj + 1 X 3 X 5 + X 6 Xs + 1 

£4 X5 xix 6 X2X7 x 3 x$ 
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These form a Poisson subalgebra with brackets { Ji, J2} = — 2 Ji J 2 , { J\, J 3 } = — 2J\ J 3 +2, with all other 
brackets following from the cyclic property and skew-symmetry. This is the algebra of Pq , after the 
permutation a : (1, 2, 3, 4, 5) 1— ► (1, 4, 2, 5, 3). It provides three commuting functions, 

lo = Ji + J2 + J3 + J4 + J5, Xi = Ji J4J2 + J2J5J3 + J3J1J4 + J4J2J5 + J5J3J1) I2 = J1J2J3J4J5, 

which (from Theorem 14. 8} are the homogeneous components of the Casimir of this subalgebra, namely 

K = trMi =I -X 1 +X 2 . 

The generators of the period 3 subalgebra are given by 

X1X3 + X 6 X S + 1 X\ + X 7 + X4,(xiX2 + x 6 x 7 ) x 2 + x s + x 5 (x 2 x 3 + X 7 X S ) 
Ai = , K 2 = , A3 = , 

x 3 x$ xiX4X 7 X2X5X8 

whose Poisson bracket relations are 

{K 1 ,K 2 } = -2K 1 K 2 +2, {K 2 ,K 3 } = -2K 2 K 3 + 2, {K x , K 3 ] = 2K X K 3 - 2. 

Up to the permutation a : (1, 2, 3) H> (1, 3, 2), this is the algebra associated with P4 , with Casimir 

K = tr Mi = -to + It, where % = K x + K 2 + K 3) X a = K 1 K 2 K 3 . 

The latter two quantities commute with each other, and we have the relation T$ — X\ + X2 + Iq — 1\ = 0. 

Since (|4.37[) holds, any four of the five functions Io,Ii,I2,XojXi provide the correct number of indepen- 
dent commuting first integrals to show Liouville integrability of the 8-dimensional map (|4.38[) . 

4.6.2 The odd case 

When N is odd, then p is odd and q is even, or vice versa. With N = 2m + 1 we find that, as for the primitives 
■Pjjm+i' the- kernel of the matrix (|4.33p is spanned by the same vector u = (1, — 1, . . . , 1, — 1, 1) T , orthogonal 
to the vectors (|4.24[) providing the symplectic coordinates yj — x Vj , as in (|4. 25[) , which are invariant under 
the one-parameter scaling group (|4.23|) . The symplectic form obtained via Lemma 12.71 gives a nondegenerate 
Poisson bracket for the yj , of the form 

{y 3 ,yk} = tjkVjVk, ±<j,k<2m. (4.39) 

In all examples we find that the Toeplitz matrix B" 1 = (ejfe) has only the entries 1, —1, above the diagonal, 
but a concise formula for these ejk in terms of the coprime integers p, q is presently unavailable. 

By Theorem I2.6[ we have an induced birational map tp in 2m dimensions, which is a Poisson map with 
respect to the nondegenerate bracket (|4.39p . We would like to assert that this is an integrable map. 

Assume for the sake of argument that p is odd and q is even. Then the J n can be written as functions of 
yj, and from ()4.39|) they should satisfy the algebra (|4.19[) corresponding to Pp+i, hence providing (p + 1)/2 
commuting integrals. The scaling-invariant quantities w n — K n K n+qi for 1 < n < q — 1, can also be written 
in terms of the yj, and (up to the permutation a) should give a Poisson subalgebra isomorphic to that of 

the functions (|4.32|) for Pg+i, providing another q/2 commuting integrals. With the constraint (|4.12j) . this 
would give m independent commuting integrals in dimension 2m, as required. 

For the rest of this section, we present examples of primitives Pjf* with odd N and q > 1. 

Example 4.16 (The primitive P 5 (2) ). For TV = 5, q = 2 the matrix (|4~33)) has null vector u = (1, -1, 1, -1, 1) T , 
and im B is spanned by the vectors 

Vi = (1, 1, 0, 0, 0) T , v 2 = (0, 1, 1, 0, 0) T , v 3 = (0, 0, 1, 1, 0) T , v 4 = (0, 0, 0, 1, 1) T . 
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Upon applying Lemma [2~71 the Toeplitz matrix B~ x is specified by its first row, namely (eij) = (0, 0, 1, 1), 
which determines the components of the nondegenerate Poisson bracket (|4.39|) for the variables yi — 
XiXi+i, i = 1, . . . ,4. Explicitly (for indices j < k) this is just 

{2/1,2/3} = 2/12/3, {2/1,2/4} = 2/12/4, {2/2,2/4} = 2/22/4, 
all other brackets being zero. The Poisson bracket is invariant under the induced map 

<P- (2/1,2/2,2/3,2/4) ^ (2/2,2/3,2/4,2/22/4(2/3 + l)/(2/i2/3))- (4.40) 
The period 3 functions J n take the form 

T Xl + 25 T X\X2 + X3X4 + 1 j X2X3 + X4X5 + 1 

J\ — , Jn — , J3 — ■ 

X 3 X1X4 x 2 x 5 

Being invariant under the scaling symmetry (|4.23[) . they can also be written in terms of the variables as 
T _ 2/12/3 + 2/22/4 T _ 2/2(2/1 +2/3 + 1) T _ 2/3(2/2+2/4 + 1) 

J\ — , J2 — , ^3 — • 



2/22/3 2/12/3 2/22/4 

The Poisson brackets between these functions follow by the cyclic property from { J l5 J 2 } = 1 — J1J2, which, 
up to rescaling by a factor of 2 and applying the permutation a : (1, 2, 3) 1— >■ (1, 3, 2), is the bracket (|4.19[) of 
the J n for P4 . The J n subalgebra provides a pair of first integrals in involution, namely 

= Ji + J2 + J3, T-\ = J1J2J3, 

which is sufficient for the map (|4.40|) to be Liouville integrable. (Since these functions are totally symmetric, 
not just cyclically symmetric, the permutation a plays no role.) 

The period 2 quantities K n , which are not invariant under the scaling (I4.23[) . are 

X\X 2 + 24X5 + I T , Xi + X 5 + X 3 (xiX 2 + XiX 5 ) 

K\ = , K 2 = . 

X2X4 X!X 3 X 5 

As for the case of P 3 in Example I4.11[ the product K\K 2 is invariant under the scaling symmetry, so can 
be written in terms of the variables y,. In fact, from (|4. 1 2[) we have K = I\ — Iq = K\K 2 — 2, where (by 
Theorem 14.81) JC is the Casimir of the bracket for the J n . 



Example 4.17 (The Case P 7 (2) ). In the case N = 7, q = 2, the Toeplitz matrix B 1 is specified by its first 
row, namely (eij) = (0, 1, 1, 0, 0, 1), which determines the components of the nondegenerate Poisson bracket 
(|4.39|) for the variables yi = XiXi+i, i = 1, . . . , 6. This bracket is preserved by the 6-dimensional map 

, , , , 2/22/6(2/32/5+2/4) 

V- (2/1, •• -,2/6) i-> (2/2, •■•,2/7), 2/7 = 



2/12/32/5 

which (by Theorem I2.6[) is induced from the map ip defined by the recurrence (|4.ip with p — 5, q = 2. 
The functions J n , which cycle with period 5 under the action of ip, can be written in terms of yi thus: 

T _ 2/12/3 + 2/22/4 j _ 2/22/4 + 2/32/5 T _ 2/32/5 + 2/42/6 

•J\ — , Jl — , ^3 — , 

2/22/3 2/32/4 2/42/5 

T 2/22/4 + 2/12/32/4 + 2/22/32/5 T 2/32/5 + 2/22/42/5 + 2/32/42/6 

J4 = , J5 = ■ 

2/12/32/5 Vi 2/4 2/6 

The Poisson subalgebra generated by these functions is specified by 

{^1,^2} = </i</2, {Ji, J3} = J1J3 — 1, 
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with all other brackets following from the cyclic property and skew-symmetry. By rescaling by a factor of 2 
and applying a permutation to order these functions as (Ji, J3, J5, J2, J4), this is seen to isomorphic to the 
algebra of the J n for Pg , so we find the commuting functions 

Xq = Jx + J3 + J5 + J2 + J47 2-1 = J1J3J5 + J2J4J1 + J3J5J2 + J4J1J3 + ^5^2>/4, 2^2 = Jl'HJz-^Ji- 

Of course, the ordering is unimportant for the totally symmetric functions Zo and 12- 
The period 2 quantities, which scale as A"i — > A 2 A"i, K 2 — > \~ 2 K 2 under (I4.23p . are 

r . X2 + X 6 + X A (xiX 2 + X 6 X 7 ) X ± X 3 + XxX 7 + X 5 X 7 + X 3 X 5 (xiX 2 + X 6 X 7 ) 

Ki = , K 2 = . 

X 2 X4Xq X\X 3 X$X 7 

From (|4.12p . the scaling-invariant combination K\K 2 can be written in terms of yi, via 

KiK 2 -2 = K.=Xo-Xi+l2- 

Example 4.18 (The Case P 7 (3) ). For N — 7, q — 3, the first row of the Toeplitz matrix B^ 1 is (e ltj ) = 
(0, 0, 0, 1, 1, 0). This defines the nondegenerate Poisson bracket (I4.39|) . which is invariant under the map 

, \ . / s 2/22/42/6(2/4 + 1) 
(yi,---,2/6) !-> (2/2, - - - , 2/7), 2/7 = , 

2/12/32/5 

in terms of the variables = XiXi + \, i — 1, ... ,6. 

The functions J n , with period 4, are not invariant under (|4.23p . but w n = J n Jn+i, n = 1, . . . ,4 are: 

(1 + 2/1 + 2/4) (2/12/32/5 + 2/22/42/e) (1 + 2/1 + 2/4) (1 + 2/2 + 2/5) 

»i = , w 2 = , 

2/12/32/42/5 2/12/5 

(1 + 2/2 + 2/s)(l + 2/3 + 2/e) W1W3 
w 3 = , w 4 = , 

2/22/6 w 2 

of which only three are independent (since 101103 = 102 104). By periodicity and skew-symmetry, all of their 
Poisson brackets follow from 

{wi,w 2 } =«>i + w 2 - w\W2, {101,103} = 104 — w 2 . 

By applying the permutation a : (1,2,3,4) H> (1,4,3,2), this is seen to be isomorphic to the algebra for 
P5 , as in Example 14.121 Hence we have two functions in involution, namely 

I\ = W\ + W4 + W 3 + W 2 , X 2 = 10i10 3 , 



where the ordering is unimportant since both functions X\ and X 2 are invariant under a. 

The necessary third function in involution is derived from the quantities K n , which cycle with period 3. 
Being invariant under the scaling (|4.23p . they can be written in terms of the variables yc. 

v 2/3(2/1 + 2/5 + 1) „ 2/4(2/2 + 2/6 + 1) 2/12/32/5(2/3 + 1) + 2/22/42/6(2/4 + 1) 
A x = , K 2 = , A 3 = . 

2/22/4 2/32/5 2/12/32/42/6 

They generate a three-dimensional Poisson subalgebra with the same relations as for the subalgebra in 
(up to a factor of 2), i.e. 

{K 1 ,K 2 } = K 1 K 2 -1, {K 2 ,K 3 } = K 2 K 3 -1, { K x , K 3 } = -K x K 3 + 1. 

There are two first integrals that are commuting functions defined on this subalgebra, which we denote by 

±0=^1+^2+ K 3 , ±i = KiK 2 K 3 , 

and the joint Casimir of the two subalgebras is given by K, = 2 — X\ + X 2 = X\ — Xq. Since {Ki, Wj} = for 
all i,j, we may use any three of the functions Xi,X 2 ,Xq,Xi to show Liouville integrability. 
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5 Linearisable recurrences from P, 



(?) 



2m 



r 2m 1 1 



(m—q) 
2(m-q) 



quivers 



In this section we consider the family of recurrences (|3.14p . as in case (iii) of Theorem 13. 121 which come from 
the quivers of the form Pj 



(</) 



j(m) 
2m 



P 



(.m—q) 



Xn+N Xn 



It is convenient to rewrite each recurrence as 
p + q = N = 2m, 



2(m-q) 
Xn-\-p Xn+q ~t~ Xn-\-mi 



(5.1) 



which (for fixed m) gives a different recurrence for each q = 1, 
by 

"7 \T 



2m 



~m — q 
T 2(m-g)> 



1. The associated matrix B is given 

(5.2) 



where tm is defined in (|4.17p . and T 2 ( m -q) denotes the N x N matrix obtained by adding q left and right 
columns and upper and lower rows of zeros to the 2(m — q) x 2(m — q) matrix r 2 ( m _,j). The simplest examples 
of the quivers corresponding to such B are shown in Figure 1. 

Observe that if gcd(m, q) = r > 1 then the quiver consists of r disjoint copies of the same type of quiver, 
but with the parameters q and m replaced by the coprime integers q/r and m/r, respectively, and similarly 
(|5.ip decouples into r copies of the corresponding recurrence. Therefore we shall assume that gcd(m, q) = 1 
from now on. With this assumption it follows from p + q = 2m that gcd(p, q) = 1 or 2 only. The case 
gcd(p, q) = 1 has a very similar structure to that of the primitives P^j ' for even N, but the case gcd(p, q) = 2 
has several new features, so we will need to distinguish between these two cases in due course. 

The family (|5.ip has some basic properties that are analogous to those of the family in case (ii) of Theorem 
13.121 as described in the previous section. In particular, all of the recurrences (|5.1j) are linearisable, and 
they have two sets of periodic functions, with periods p, q respectively, which lead to the construction of first 
integrals. 

5.1 More linear relations with periodic coefficients 

By analogy with (|4.2p , the recurrence (|5 . 1 [) can be written in the form 



det 



(5.3) 



•En •E-n+q 
^n+p %n-\-N 

The above identity is the relation for a 2-frieze [3S] , and it implies that the iterates of (|5.ip can be placed in 
the form an infinite 2-frieze. Using Dodgson condensation once again to condense a 3 x 3 determinant, we 
have 

det ^ n — {^n+m^n+N -\-m ^n+m+g^n+m+p) I '-^n+N 1? 

with \l/„ as in (14. 3p . Then condensing the appropriate 4x4 matrix A„ in terms of 3 x 3 minors yields 



det A r , 



x r , 

%n+2p 
X 



•Ejl+q 
Xn+N 
Xn+N+p 



•Kn-\-2q %n+3q 
Xn+N+q Xn+N+2q 
X n +2N X n +2N+q 



n+3p X n +N+2p X n +2N+p X n +3N 



0. 



As in the case of the primitives considered before, the left and right kernels of the singular matrix A„ yield 
linear relations between the x n - 



Lemma 5.1. The iterates of the recurrence J[ ) satisfy the linear relations 

Xn-\-3q Jn+m X n +2q ~t~ Jn Xn+q X n 0, 
Xn-\-3p Kn+m Xn+2p K n Xn+p Xn — 0, 

whose coefficients are periodic functions of period p, q respectively, that is 

Jn+p = Jn, K n +q = K n , f ' OT all n. 



(5.4) 
(5.5) 
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Proof: Upon solving A n k n = 0, using the first three rows, it is convenient to normalise the first entry of k„ 
to be — f , and solve a 3 x 3 system to find the other three entries. Then from Cramer's rule and det ^ n = 1 
the fourth entry must be +1, so that this vector has the form k„ = (— 1, J„, J„, 1) T . The 3x3 linear system 
coming from the first three rows of the equation A„ +p k n+p = is the same as that coming from the last 
three rows of A„k„ = 0, which implies that J n + P = J n and J ra + P = J„. Applying Cramer's rule in the first 
three rows of A n k„ = together with Dodgson condensation also implies that 



and similarly 



X n +2r. 



Xn X n J r 2q •Kn J r3q 

X n +p x n+N+q x n+N+2q 
X n +2p %n+2N Xn+2N+q 



En+ 2M *£n+m+p 



l 



where E n = 



•&n •£n-\-2q 
%n-\-p %n-\-N-\-q 



(5.6) 



(5.7) 



Permuting the first and second columns of the determinant in the identity det ^ n = 1 and expanding in 
terms of 2 x 2 minors, and then doing the same thing after permuting the second and third columns instead, 
leads to the formulae 



Ffi X rL J rTn 



E 



n-\-p ^"n-\-m-\-p 



(5.8) 



respectively. Then the combination x n +2m x n+5m-p(Jn+m + Jn) can be rewritten as a sum of determinants, 
whose entries can be expanded using each of the identities (|5.8p , to yield 



^n+2m ^n+2m+2g 
^n+2m+p •£n+2m+iV+q 



n+2m 



o. 



as required. This proves (|5.4|) . and the relation (|5.5j) follows by symmetry, from considering the left kernel 
of A„. □ 



Remark 5.2. The four-term linear relations (|5.4j) and f|5 . 5|) . together with det^n = 1, should be compared 
with those of the pentagram map |31j . but there the coefficients of the second and third terms are independent. 

Remark 5.3. When q — 1 the coefficient K n has period 1, so K n+1 = K n = JC for all n, and the recurrence 
(|5.5[) is just the constant coefficient, linear difference equation 



X n +3N-3 - K. X n +2N-2 + K. Xn+N-1 ~ X n = 0. (5.9) 

An immediate consequence of the latter relation is an inhomogeneous version of (|4.6[) , namely 

x n +2N-2 - (/C - 1) Xn+N-i + x n = F n , where F n+N -i = F n , 

for some quantity F n (which has period p — N — 1 in this case). Thus, by a minor modification of Corollary 
14.51 one can find explicit formulae for x n in terms of Chebyshev polynomials. 
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Remark 5.4. When q — 2, the coefficients in the recurrence (|5.5I) have period 2, so if n +2 = ^Gi for all n, 
and (since gcd(m,q) = 1 implies that m is odd) we have K n+rn = K n+ \ 1 whence 

x n +3 P - K n+ i x n+2p + K n x n+p - x n = 0. (5.10) 

This is a four-term linear relation, whose coefficients alternate with the parity of the index n. 

The existence of the periodic quantities J n and K n means that, as for the case of primitives considered 
in the last section, one can construct first integrals by taking cyclically symmetric functions of each of these 
sets of quantities. When gcd(p, q) = 1 we find one relation between these two sets of quantities, and when 
gcd(p, q) = 2 we find two relations, which will be discussed below. 



5.2 Monodromy and linear relations with constant coefficients 



This subsection follows closely the discussion of subsection 14.21 for the case of primitives. However, it will 
subsequently be necessary to refine the discussion further, depending on whether gcd(p, q) = 1 or 2. 
The relation (I5.4[) implies that the matrix \P n satisfies 




* ?l+g = * ?l L„, L n = | I ii -J„ | . (5.11) 
On the other hand, the recurrence (|5.5p yields 

/ 1 

<F n+p - t n L„ = 1 | . (5.12) 

y 1 K'n K-n+m 

As before, upon taking the ordered product of the L ra over p steps, shifting by q each time, we have the 
monodromy matrix 

M„ := L„L„ + q . . . L„ + ( p _!) 9 = ^ n ^ n+pq . (5.13) 
Taking the ordered product of the L„ over q steps, shifting by p each time, gives another monodromy matrix 

M„ := Vn+p^n 1 . (5.14) 

From the cyclic property of the trace it follows that 

K„:=trM„ = trM n . (5.15) 

The periodicity of L„, together with (|5.13[) . implies that /C„+ p = /C„, and similarly, from (|5.14[) , we have 
K-n+q = K-n- If the periods p and q are coprime, then K, n = K, ^constant, for all n, hence K. is a first integral 
for the map if corresponding to (|5.1[) . However, if gcd(p, q) = 2 holds instead then we have JC n +2 — K-m so 
this quantity has period 2. 

Once again the general result of Lemma fe.lOl can be applied here, to show that the iterates of (|5.1[) satisfy 
a linear relation with constant coefficients. 

Proposition 5.5. The iterates of the nonlinear recurrence \5. 1)) satisfy a linear relation of order 3pq with 
constant coefficients. It is a four-term relation if gcd(j>,q) = 1, and a seven-term relation if gcd(p,q) = 2. 

Proof: From (|5.4p . Lemma 13.101 implies that x n satisfies a linear recurrence of order 3pq with constant 
coefficients, for which the gaps between indices of adjacent terms with nonzero coefficients are of size p. On 
the other hand, (|5.5[) implies a linear recurrence of the same order with gaps of size q. Thus the actual size 
of the gaps must be the lowest common multiple of p and q. Hence, when p and q are coprime, the gaps 
are of size pq, giving a four-term relation, while gcd(p, q) = 2 gives a seven-term relation with gaps of size 
pq/2. ' ' □ 

Below we provide refined versions of the preceding result, with more precise details of the coefficients, by 
considering the cases gcd(p, q) — 1 and gcd(p, q) — 2 separately. 
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5.3 The case gcd(p, q) = 1 

The discussion of the case where p and q are coprime is almost identical to that for the primitives pffl 
with N even, as in subsection 4.7 above. The integers p and q are both odd, and for the matrix (|5.2[) in 
each case we find that det-B = 4 whenever p or q is divisible by 3, and detB = 1 otherwise. From such a 
nondegenerate matrix B we get an invariant Poisson bracket of the form (|2.6p . which is specified uniquely 
(up to scale) by the Toeplitz matrix C = B^ 1 . However, in this case a general closed-form expression for the 
entries of C, analogous to the formula (|4.36[) for the Poisson bracket of the even primitives, is not available 
to us at present. 

We now consider the associated functions J n and K n , of periods p and q respectively. From the expression 
(I5.13|) . the first integral JC defined by (|5 . 1 5[) is a cyclically symmetric polynomial in the J n , n = 1, . . . ,p, and 
from (|5.14[) it is also a cyclically symmetric polynomial in the K n , n — 1, . . . , q. Thus the equality of the 
traces in (j5 . 15[) provides a single functional relation between these two sets of functions. Note that we also 
have the same phenomenon as in (|4.34p with regard to the different expressions for JC as a function of the 
quantities J„, when we compare the cases with q = 1 and q > 1, for the same value of p: up to the action of 
a suitable permutation a of (1, . . . , p), the two expressions are identical; and the analogous statement applies 
to JC considered as a function of the K n . 

Now observe that all of the preceding comments concerning JC apply equally well to the quantity 

JC := trM; 1 = trM" 1 , 

which is also a first integral (so this definition holds for any n). We would like to assert that in fact JC = JC. 
To see this, note that trL„ = K n+m , and trL" 1 = K n . Thus, in the case q=l, when K n = JC = constant, 
we have JC = trM„ = trM,^ 1 = JC by (|5.14[) . This then implies that, in terms of functions of Ji, . . . , J p , we 
have 

trM„ = trM" 1 (5.16) 

for q — 1, and clearly this identity remains true when q > 1, since (for fixed p) the functions JC and JC are 
obtained from the case q = 1 by applying the same permutation a to both sides. This allows us to make a 
more precise statement than Proposition 15.51 



Theorem 5.6. When gcd(p,q) = 1, the iterates of the nonlinear recurrence H5.1\) satisfy the linear relation 

•En+Zpq JC X n J r 2pq ~\~ JC Xn-\-pq — Oi (^'-^^) 

where JC is the first integral defined by \5.15\) . 

Proof: Using (|5.13|) we see that ^ n + pq = *nM„, so ^ n +2 Pq = ^ n + P qM n M n+pq = *„M^, by periodicity 
of M„, and similarly ^„ +3p( j = ^M^. Applying Cayley- Hamilton to both M„ and M" 1 , and noting that 
detM„ = 1 and trM„ = JC, as well as f|5 . 16[) . yields 

Vn+3p q - JC^ n+2pq + JC^ n + pq - = *„(M3 - JCMl + JCM n - I) = 0. 

The (1,1) component of this equation is just (I5.17[) . □ 

Remark 5.7 (Another construction of the integral for P^). We may think of the map ip in Example 14. Ill as 
coming from the recurrence y n +iVn — Un+i + Un+i, which is of the same form as (|5.1[) , with p = q = m = 1, 
although not directly obtained from a cluster mutation. Replacing x n — > y n and JC — > fC in (|5.17[) . and 
solving, leads to the first integral 

£ = Vn+3 ~ Vn 
Vn+2 — Un+1 

When this is written purely in terms of y n ,y n +i, using the map, we find the previously obtained quantity 

T , , (y n + y n +i)(y n + y n +i + l) 

VnVn+l 
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For the Liouville integrability of the map <p corresponding to (I5.1[) . the counting of first integrals appears 
to be the same as for the even primitives P^m- The two sets of quantities J„ and K„ generate Poisson 
subalgebras of dimensions p and q, which should contain [p + l)/2 and (q + l)/2 commuting integrals, 
respectively, including the Casimir K, which is common to both subalgebras. Taking the constraint (I5.15[) 
into account, this produces m = (p + q)/2 independent commuting integrals in terms of the Xj, as required. 

However, while the above counting argument is plausible, it rests on some unproven assumptions. Every- 
thing relies on the structure of the Poisson bracket for the J n in the case q = 1, since all other subalgebras 
of J n or K n should be isomorphic to one with q = 1. Yet in general (excluding m = 2), this Poisson bracket 
consists of three homogeneous parts (of degrees 0, 1 and 2), which (unlike the bracket (14. 19)) for the prim- 
itive -P^m) does not have an obvious splitting into a bi-Hamiltonian pair. Moreover, while every cyclically 
symmetric polynomial function of the J n is a first integral of the map ip, we do not yet have an algorithm 
for selecting an involutive set of them. 

Remark 5.8. In [21] a quadratic Poisson structure is presented for the coefficients of the four-term linear 
relations for twisted polygons, together with corresponding monodromy matrices and commuting first inte- 
grals for the pentagram map. However, a Dirac reduction of this bracket to the case of (|5.4I) or (I5.5[) gives 
only the trivial bracket. A general approach to Poisson structures related to twisted polygons is described 
in [27], which should shed some light on the situation here. 

For want of more general statements, we illustrate the foregoing discussion with several examples of the 
intcgrable systems that arise in this case. 

Example 5.9 (The quiver P 4 (1) - P 4 (2) + P 2 (1) ). This example was first studied in an ad hoc way in |21) . 
The appropriate matrix (|5.2I) . which corresponds to the first quiver in Figure 1, is obtained by setting c = 1 
in (|2 . 14[) . The explicit form of the map is 

(p : (x 1 ,x 2 ,x 3 ,x i ) i-> (x 2 ,x 3 ,x 4i ,x 5 ), x 5 = X2Xi + X3 . (5.18) 

Xi 

The invariant Poisson bracket (I2.6[) for this map is given by the Toeplitz matrix C = 2B^ 1 , with top row 
(ci J) = (0, 1, 1, 2). The period 3 functions Ji take the form 

x 3 (x 2 + xix 3 ) + x±(x\ + x\) _xix i + x\+x\ x 3 (x 2 + X1X3) + X4(xiX4 + x%) 
J\ — , J2 — , ^3 — • 

X\X 2 Xi X 2 X 3 X1X3X4 

The Poisson brackets between these three functions follow by the cyclic property from 

{Ji,J 2 } = JiJ 2 -2J 3 . (5.19) 
This example is exceptional, in that the bracket for the Ji is the sum of only two homogeneous terms: 



P = P (2) + p(i) ; where p(2) I _ J± j 2 J 2 J 3 , pW = 





Each of the tensors specified by and P^ 2 ) satisfies the Jacobi identity, so (since their sum is a Poisson 
tensor) they define a compatible pair of Poisson brackets. 
The first integrals 

Hi = Ji + J 2 + J3 , H 2 = J\J 2 J 3 

satisfy the bi-Hamiltonian ladder 

PiVfti = 0, P2VH1 = P1VH2, P2VH2 = 0, 
so they commute with respect to the bracket defined by (|5.19|) . The quantity 

K = 3-Hi+H 2 (5.20) 
provides the Casimir of this bracket; following (|4.34j) . we will find it useful to denote this quantity by £3,1. 
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Example 5.10 (The quiver P 6 (1) - P 6 (3) + P 4 (2) ). Mutation of the second quiver in Figure 1 gives the map 

X 2 X 6 + X 4 

if. (Xi,X 2 ,X 3 ,X4,X 5 ,X 6 ) H> (x 2 ,X 3 ,X 4 ,X 5 ,X 6 ,X 7 ), x 7 = . 

Xi 

The invariant Poisson bracket for this map is given by (|2.6|) . with the coefficients specified by the Toeplitz 
matrix C = B^ 1 with top row {c\j) = (0, 0, f , 0, 1 , 1). The functions Ji, which cycle with period 5 under the 
action of tp, take the form 

xxx 5 + x 2 x e + x 3 x 4 _ x 4 (x 3 + XxXs) + (xi + x 3 )x 2 x 6 x 3 (x 2 + x 4 ) + xix 5 ) 
Jl — 1 <H — ) ^3 — > 

X 2 X 5 XiX 3 X 6 X2X4 

x 4 (x 3 + x 5 ) + x 2 x 6 ) x 4 (x 3 + X1X5) + (xix 5 + x 2 x 3 )x 6 

•Ji — , Jb — ■ 

X 3 X5 X1X4X6 

The Poisson bracket between these five functions follow by the cyclic property from 

{Ji, J2} = -J1J2 -Ji + l, { Ji, J 3 } = 2Ji J 3 . 
This Poisson bracket is the sum of three homogeneous terms, 

P = PW+PW+P(0), with PL 2) = 4 2) J, J fe , P^=c«J, +2 , P^=c^, 
where c ik are the Toeplitz matrices with top rows given by 

(<S) = (0,-1,2,-2,1), (c^) = (0,-1,0,0,1), (eg) = (0,1,0,0,-1). 

Both P^ ^ and P^ 2 ^ satisfy the Jacobi identity, but P^ 1 -* does not, so we cannot think of this sum as some 
sort of Poisson compatibility. 

The Casimir for the 5-dimensional Poisson algebra generated by the Ji is the trace of the monodromy 
matrix, as in (|5.15j) . It can be written as the sum 

IC = n 1 -H 2 +H 3 , (5.21) 

where each of the components is a first integral: 

55 5 

=~^^{Ji — JiJi+i), %2 = / Jj+lJj+2 ~ JjJj+lJi+2), H 3 — 

i—1 i—1 i—1 

We find that {Hi,Hj} = for all so the 6-dimensional Poisson map <p has the correct number of first 
integrals in involution. For later comparison, we denote K. in (|5.2ip by K.^ t \. 

Example 5.11 (The quiver P^ 1 ' — + Pg 3 '). The map obtained from mutation of this quiver is 

x 2 x 8 + x 5 



tp: (xi,X2,x 3 ,X4,X5,x 6 ,x 7 ,x 8 ) n- (x 2 , x 3 , x 4 , x 5 , x 6 , x 7 , x 8 , x 9 ), x 9 



The corresponding non-singular matrix B defines an invariant Poisson bracket (12 .6[) with matrix C = B^ 1 . 
The top row of the Toeplitz matrix C is (c.\ j) = (0, 1, 0, 0, 1, 1, 0, 1). The functions Jj with period 7 can be 
determined from Ji, which takes the form 

XiX 6 + X 2 X 7 + X 3 X 5 

J\ — ! 

X 2 X6 

The remaining six functions are obtained by applying ip* Ji — Ji+i, with {tp*) 7 Ji = Jj. The Poisson bracket 
relations between these functions follow by the cyclic property from 

{Ji, J2} = 2Ji J2 — J5, { J\, Js} = — J\J 3 + 1, {Ji, J4} = — J\ Ji- 
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Again, this Poisson bracket is the sum of three homogeneous terms, P = P^ 2 - 1 + P^ + P(°), where p(°) and 
P^ 2 ) satisfy the Jacobi identity, but P^ does not. 

The Casimir of the Poisson subalgebra generated by the Ji is again K. (the trace of the monodromy 
matrix), but in this case it is not clear how to split the Casimir into four pieces that Poisson commute, 
as required for Liouville integrability. Of course, there are still seven functionally independent invariant 
functions (built from cyclically symmetric combinations of the J,) and we expect that four commuting 
functions exist. 




4 3 

(a) Example 15.91 




5 " 4 

(b) Example IBTTOl 

Figure 1: The first three quivers in this class. 




5 4 

(c) Example 15. 161 



Example 5.12 (The quiver P^ — P 8 W + P^ 1 )- For this quiver, with p = 5 and q = 3, the recurrence is 

The corresponding birational map tp is Poisson with respect to a log-canonical bracket (|2.6|) defined by the 
inverse of the associated matrix B. Upon setting C = 2i3~ 1 , the bracket is given by the first row of this 
Toeplitz matrix: (c 1}j ) = (0, 1, 1, 0,-1, 1, 2, 1). 

Since gcd(p, q) = 1, we have the same phenomenon as in the discussion of (|4.34[) . From (15.13[) . we have 

M n := L„L n _|_3L n+ 6L„ + 9L n+ i2 = L„L„ + 3L„ + iL n _|_4L„ + 2, 

since L„ has period 5, so the permutation a : (1, 2, 3, 4, 5) i-> (1, 4, 2, 5, 3) gives the trace as 

= f^5,3[Ji, J2, J3, J4, J 5) — fc-5,i(Jii J4, Ji, J5, J3), 
where K,$ t i is given by the formula (|5.21j) in Example 15.101 Similarly, from (|5 . 14[) we find 

M„ := L„ + ioL n+ 5L„ = L„ + iL„ + 2L n , 
where L„ has period 3, so with the same notation as in (|4.35j) we also have 



>(4) 



,(1) 



where K. 3i ± is given by (|5.20l) in Example l5.9l Since £3,1 is totally symmetric, the permutation a : (1, 2, 3) 1— > 
(1, 3, 2) makes no difference to the result. 

The explicit forms of the period 5 functions Ji, appearing as entries in the matrices L„, all derive from 



Jl 



£1X3X8 + X3X5X7 + X4'X 6 X 8 + 24X7 

232428 
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by acting with the map ip. From the above bracket for the Xj, determined by the matrix C, the Poisson 
brackets between the Ji can be calculated as 

{Ji, J 2 } = 4Ji J 2 , { Ji, J 3 } = 2Ji J 3 + 2J 2 - 2, 

with all other brackets being deduced through the cyclic property. The resulting 5-dimensional Poisson 
subalgebra is isomorphic to that for the J, in Example 15.101 as can be seen by applying the permutation a 
and rescaling by an overall factor of 2 (which depends on the choice of scale for c^). Therefore, subject to 
this permutation, it follows that the same three functions Hi are in involution: 

5 5 5 

"Hi = — JiJi+3), H 2 = / l(JjJj+lJi+3 ~ JiJi+lJi+z)) % = TT Ji- 

i—1 i—1 i—1 

The required fourth function that commutes with these Hi is derived from the algebra of the Ki, which 
appear as entries in the matrices L„. In terms of Xj, we have 

Xi x 6 x 2 x 5 x$ x 2 x s 
ill = 1 1 1 1 1 , 

£6 X\ £3X6 X\X± X4X7 

with K% = ip* K\ and K% — {ip*) 2 K\ providing the other two functions which cycle with period 3 under the 
action of the map. The Ki generate a 3-dimensional Poisson subalgebra, whose brackets are all determined 
by acting with <p on the single relation 

{K 1 ,K 2 } = -K 1 K 2 + 2K 3 . 

This algebra is isomorphic to that of the Jj in Example 15.91 (as is seen immediately by applying the permu- 
tation a), so it contains the two commuting quantities 

Hi = K\ + K\ + K\, H 2 = K1K2K3. 

The subalgebras of Ji and Ki share the joint Casimir JC = Hi — 1L 2 + V.3 = 3 — Hi + H2, which gives a 
single relation between these two sets of functions. Since {Ki, Jj} = 0, for all i, j, we can take any four of 
the first integrals ^1,^2,^3,^1,^2 as a commuting set, which demonstrates that this 8-dimensional map 
ip is integrable in the Liouville sense. 



5.4 The case gcd(p, q) = 2 

The discussion of the case where p and q are both even involves some new features. Upon setting p = 2p, 
q = 2q with gcd(p, q) — 1, the fact that gcd(m, q) = 1 implies that m = p + q is odd, hence either q is 
odd and p is even, or vice versa. For the matrix (|5.2[) in each case we find that det B = whenever p or 
q is divisible by 3, and det B = 9 otherwise. These two possibilities lead to quite different behaviour, so 
eventually we shall have to distinguish between them. For the time being we concentrate on the associated 
linear relations, which do not depend on whether B is degenerate or not. 

Since q is even, observe that the ordered product of matrices L„ in (|5.13l) now cycles only through indices 
with the same parity, so the product cycles twice through p = p/2 terms. Thus we see that M„ is a perfect 
square, and it is convenient to take the square root 

My 2 = L„L n+9 . . . L ri +(p_ i) g , 

and similarly for M„ . We know that in this situation the quantity K. n in (|5.15[) has period 2, but for our 
purposes it is more useful to consider the quantity 

/C;:=trMy 2 =trMy 2 , (5.22) 
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which cycles with period 2 for the same reasons. 

In this setting, with gcd(p, q) = 2, the algebraic structure in terms of the functions J„ and K n is based 
on that for the case q — 2 (up to suitable permutations), similarly to the way that for gcd(p, q) = 1, and for 
the primitives, this structure is based on the case q = 1. 

Now when q = 2 we have the quantities K n with period 2, giving 

tr My 2 = tr L n = K n+1 , tr M^ 2 = tr 1 = if n) 

and hence the identity 

trM j ; 1 /2 =trMj 1 / + 2 1 =/C; +1 (5.23) 

1/2 

holds. For even q > 2, with p fixed, the formula for trM„ as a function of J\, . . ., J p is identical to that 
for the case q = 2, up to a permutation, which means that the formula (|5.23p holds in general. 

Theorem 5.13. When gcd(p,q) — 2, the iterates of the nonlinear recurrence i5.1\) satisfy the linear relation 

x n+3pq/2 ~ K-n x n+pq + K-n+1 x n+pq/2 — x n = 0, (5.24) 

where /C* is the period 2 quantity defined by h5.2'2\) . 



Proof: This follows by essentially the same argument as in the proof of Theorem 15.61 applying Cayley- 
Hamilton to M^ 2 and M^ 1 ^ 2 , and making use of (|5.23[) . □ 

Remark 5.14. With respect to the functions J„ and K n , the main new feature in this case, compared 
with the case gcd(p,g) = 1, is that here we have two quantities JC\ and K,\, so the identity (|5.22l) for 
n = 1,2 provides two independent relations between these two sets of functions. Given that these are the 
only relations, the cyclically symmetric functions of the J„, and those of the K n , together provide N — 2 
independent first integrals. 



We are now ready to present a further refinement of Proposition 15.51 

Theorem 5.15. When gcd(p,q) = 2, the iterates of the nonlinear recurrence H5.1]) satisfy the constant 
coefficient, linear relation 

x n+3pq ~ 13 x n+5pq/2 + C X n -\-2pq — T> X n +3pq/2 + C X n + pq — B X n +pq/2 + x n = 0, (5.25) 

where the non-trivial coefficients are first integrals, given by 

B = Kl+K* 2 , C = KIK* 2 + Kt + K* 2 , V = (/q) 2 + {JC* 2 f + 2. 

Proof: For the sake of argument, suppose that q is odd and p is even. To get the seven-term relation (|5.25p . 
we consider the 6x6 matrix that is specified in terms of its (j,k) entry by $„ = { x n +pq(j-i)+q{k-i))- The 
linear recurrence (|5.4p implies that 



$ n+I j = $„L*, where 



/ 1 \ 

1 

1 -J n 

1 

1 J n+m 

\ 1 



The monodromy matrix corresponding to p iterations of the latter equation is 

IVf* T * T * T * 

lvl n — ^n+q • • • ■ Ll n+(p-l)(jJ 
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so that 

$„+^ = $„M;. (5.26) 
We need to show that the characteristic polynomial of the 6x6 matrix M* is given by 

X (Q = ( 6 -B( 5 +C( 4 ~V( 3 +C( 2 -B( + l, (5.27) 

as once we have determined this the recurrence (I5.25P follows immediately, by applying the Cayley-Hamilton 
theorem to (|5.26[) . 

It is convenient to conjugate all of the 6x6 matrices by the permutation matrix corresponding to 
(1, 2, 3, 4, 5, 6) — > (1, 3, 5, 2, 4, 6), which reduces everything to calculations with 3x3 blocks. This gives 



1 / ' ^n^n+q" I o L 



n+q 



so that that for the monodromy matrix we have 



L n+ qL n+ g +q . . . L n+q+ (p_ 1 - )q J \ o 



1/2 

IV /T 

n+1 



1/2 



Thus the characteristic polynomial of M* factors as the product of the characteristic polynomials of M 
and M^, 

x(0 = (C 3 - icic 2 + K.1 C - i) (C 3 -ic* 2 c 2 + icic- i), 

which multiplies out to give (|5.27[) with the correct coefficients B,C,T>. 

An analogous argument holds when q is even and p is odd. □ 

To discuss the Liouville integrability of the systems that appear when gcd(p,q) — 2, it is necessary to 
give a separate treatment according to whether the matrix B is invertible or not. 



5.4.1 Nondegenerate B matrix 

When matrix B is nondegenerate, it appears that the Liouville integrability of the symplectic map ip should 
follow by very similar arguments to those for the case gcd(p, q) — 1 (or for the primitives with even iV). The 
map preserves a Poisson bracket of the form (|2.6p , which is specified uniquely (up to scale) by the Tocplitz 
matrix C = B~ x . There are the two sets of quantities, J„ and K n , with periods p and q respectively, and by 
Remark 15.141 these provide N — 2 independent first integrals, but it is necessary to find m = N/2 integrals 
in involution. 

Given that the J„ generate a Poisson subalgebra of dimension p, and that both of the quantities JC{ and 
K,2 defined by (I5.22p are Casimirs, with the symplectic leaves being of dimension p — 2 = 2(p— 1), a further 
p—1 independent commuting functions of the J n are required in order to define an integrable system on this 
subalgebra. Similarly, given that {Ji,Kj} = for all and the K n produce a g-dimensional subalgebra, 
with the same functions K.\ and K,2 as Casimirs, it is necessary to have an additional q — 1 independent 
functions that define an integrable system in terms of K n alone. The quantities B and C which appear as 
the coefficients in the linear relation (|5.25|) are first integrals, as well as being joint Casimirs for the two 
subalgebras, so combining these with the two sets of first integrals gives a total of 2 + (p — 1) + (q — 1) = m 
independent functions in involution, as required. 

Since we do not have a general proof of all the foregoing assertions, here we will merely illustrate the 
discussion with the simplest example of this kind, which arises for p — 4, q = 2. 

Example 5.16 (The quiver P 6 (2) - P 6 (3) + P 2 (1) ). Mutation of the third quiver in Figure 1 gives the map 

ip: (x l ,x 2 ,x 3 ,x 4: ,x 5> x 6 ) k> (x 2 ,x 3 , x 4 ,x 5 ,x 6 ,x 7 ), x 7 = X3X5 + Xi . (5.28) 

Xi 
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The inverse of the non-singular matrix B defines a Poisson bracket (|2.6I) , invariant with respect to ip. Taking 
C = SB^ 1 , the top row of this Toeplitz matrix is (cij) = (0, 1, 1, 0, 2, 2). 

The period 4 functions Jj, generated from Ji by the action of the map, are 

_ XlX 2 X 6 + X2X4X5 + £3X4X6 -(- X3X5 _ XlX 6 + X 2 X 3 + X4X5 
Jl — ) J2 — ! 

X 2 X 3 X 6 X3X4 

X1X3X4 + X2X3X5 + xix^Xq X1X3X5 + x 2 x 4 x 6 + xxx 2 x±x§ + xix 3 x 4 x 6 + X2X3X5X6 

J 3 = , J4 = . 

X1X4X5 X1X2X5X6 

The Poisson brackets between these functions all follow by the cyclic property from 

{ Ji, J 2 } = 3Ji J 2 - 3, { J u J 3 } = 3( J 2 - Ji). (5.29) 

From the first equality in (|5.22|) we have the two functions 

/C 1 = J 2 J4 — Ji — J3, /C 2 = Ji J3 — J 2 — J4, 

which are both Casimir functions of the 4-dimensional algebra generated by the J i; but cycle with period 2 
under the map. Symmetric functions of K,\ and 1C 2 are first integrals as well as being Casimirs, so we may 
take the quantities 

B = jq + K%, C = K*K 2 + B, 

as in Theorem 15. 151 The homogeneous components of B are automatically first integrals, and they Poisson 
commute: 

B = -n 1 + n 2 

where 

Hi = Ji + J 2 + J 3 + Ji, H 2 = Ji J 3 + J 2 Ji, with {Hi,H 2 } = 0. 

Hence, we may take %i,£>,C as three independent first integrals in involution, which proves Liouville in- 
tegrability of the 6-dimensional map (|5.28[) . Clearly there are other choices, for instance Hi,H 2 ,C, which 
would do just as well. 

Remark 5.17. Note that in the preceding example, since q = 2 we have K\ — JC 2 and K 2 = JC^, which 
Poisson commute with each other, so the subalgebra of the Ki is trivial. 

5.4.2 Degenerate B matrix 

When B is degenerate, which only happens when either p or q is a multiple of 6, the maps that arise have 
features that make them much more like the odd primitives P 2 m+i than the other cases with even N. 

For these particular cases, the matrix B has a two-dimensional kernel, which is spanned by two integer 
vectors of the form 

m = (1, 1, 0, -1, -1, 0, . . .) T , u 2 = (0, 1, 1, 0, -1, -1, . . .) T , 

where in each vector the components continue to repeat the same blocks of six numbers, until the final block 
which is truncated (of length 2 or 4, since TV = 2m is even and not a multiple of 6). Hence there is a 
two-parameter scaling group, which acts by 

(xi,X2,x 3 ,X4,x 5 ,x 6 , . . .) -> (Axi,A^X2,AiX3,A _1 x 4 ,A _ V _1:;c 5,M _1 a ; 6,...), (A,m) € (C*) 2 . (5.30) 

This action extends to all the iterates x„ of (|5.1j) ; the pattern repeats itself on each successive block of six 
adjacent iterates. 

Now im B is spanned by 

v j ■= e j ; - e j+i + e j+2, j = 1, ...,2m -2, (5.31) 
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where Bj is the jth standard basis vector. Hence, by Lemma \2. 71 the coordinates 

% = ^, j = l,...,2m-2, (5.32) 

are invariant under the scaling (|5.30l) . and the degenerate form (|2.10p pushes forward to a symplectic form 
(|2.22|) in dimension 2m — 2. The coefficients of the latter are obtained from a skew-symmetric matrix 
B = (bjk), whose inverse provides a nondegenerate Poisson bracket for the yj, i.e. 

{yjiVk} = ejkVjVk, l<j,k<2m-2, (5.33) 

where B^ 1 = {tjk) must be a Toeplitz matrix, because the coordinates (|5.32[) transform as ip*yj — 2/j+i 
under the map corresponding to (|5.1|) . Upon applying the rest of Theorem 12.61 we see that this induces a 
symplectic map (p on the variables yj. 

To prove the Liouville integrability of all of the maps ip that arise in this way, we require a general 
expression for the coefficients ejk in (|5.33[) . which is presently lacking. However, it is possible to give a 
plausible argument for the counting of first integrals, which agrees with all examples we have checked so far. 

Suppose, for the sake of argument, that 6\q, hence 6/\p. The scaling action of (C*) 2 on the Xj, as in 
(|5.30j) . extends to an action on the coefficients J n that appear in the linear equation (|5.4j) . However, the 
indices of the terms x n +j q for j — 0, 1, 2, 3 differ by multiples of 6, which means that they all scale the same 
way, and hence the period p quantities J n are invariant under this scaling, and can be expressed in terms of 
the yj given by (|5.32l) . Given that the J n generate a p-dimensional Poisson subalgebra with respect to the 
bracket (15.33|) . with two Casimirs given by the quantities K,\ and K,\ as in (|5.22p . a further p — 1 independent 
commuting functions of the J n are needed to have an integrable system defined on this subalgebra. Similarly, 
the scaling action (|5.30|) extends to an action on the coefficients K n in (15. 5| . but now the terms x n +j P for 
j = 0, 1,2,3 do not all differ by multiples of 6, so they scale differently. This means that there is a non- 
trivial scaling action of the two-parameter group (C*) 2 on the quantities K n , n = 1, . . . , q, for which there 
should be q — 2 invariant monomials. The invariant monomial functions of K n , which we denote by Wi for 
i = 1, . . . ,q — 2, can also be written in terms of the original variables Xj. The fact that they are invariant 
under (fOP]) means that these u>i can be written as functions of the symplectic coordinates yj as well. Given 
that the Wi generate a (q — 2)-dimcnsional Poisson subalgebra in the (2m — 2)-dimensional space with the 
bracket (|5.33p . and that the quantities K\ and K,2 are Casimirs for this subalgebra too, an integrable system 
is defined on the (q — 4)-dimcnsional symplectic leaves by an additional set of q — 2 commuting functions of 
the Wi. Supposing further that {Ji,Wj} = for all i,j, we take the quantities B and C from (|5.25|) . which 
are both first integrals and joint Casimirs for the two subalgebras, and combining these with the two sets of 
first integrals gives a total of 2 + (fi — 1) + (q — 2) = m — 1 independent functions in involution, as required 
for the Liouville integrability of the symplectic map tp in dimension 2m — 2. 

Since we do not have a general proof of all the preceding assertions, here we will only present the simplest 
example of this kind, which arises for p = 6, q = 4. 

Example 5.18 (The quiver P$ — pffi + P% ). The recurrence arising from mutation of this quiver is 

X n +W x n = Xn+6 X n +4 + X n +5- (5.34) 

The matrix B has a two-dimensional kernel, spanned by the integer vectors 

Ul = (1,1,0,-1, -1,0,1,1,0, -1) T , u 2 = (0,1, 1,0,-1,-1, 0,1, 1,0) T , 

which generate the action of a two-parameter scaling group on the iterates of (|5.34p . as given in (|5.30p . 
Taking the scaling-invariant variables yj, given by (|5.32[) for j — 1, . . . , 8, we apply Lemma 12.71 to find the 
symplectic form Qj expressed in these coordinates, which leads to the nondegenerate Poisson bracket specified 
by 

{2/1,2/5} = 2/12/5, {2/1,2/e} = -2/12/6, {2/1,2/7} = 2/12/7, (5-35) 
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where all other brackets are either zero or follow from the Toeplitz property/skew-symmetry. From Theorem 
12.61 this Poisson bracket is preserved by the induced map 

I \ ,t \ 2/42/52/6(2/5 + 1) , K „„. 

v- [yi,---,ys) !->■ (2/2, •■•,2/9), 2/9 = ■ (5.36) 

2/12/22/8 

Note that (in contrast to the foregoing discussion) in this example q is not a multiple of 6, but rather p 
is. The period 6 functions J n are given by the formula 

X n Xjij-7 X n — \X n j-% 

J n = 1 1 , n = l,...,6, 

X n +4 X n +3 X n +3X n +4, 

where each of the above quantities can be written as a function of x\, . . . ,x\o by iterating the recurrence 
(I5.34[) either forwards or backwards. Since the gaps between the indices of the terms Xj in the linear relation 
(I5.4[) are not all multiples of 6, these J n are not invariant under the scaling action (|5.30[) . but instead they 
transform as follows: 

Ji — > A 2 /iJi, J2 — > X/j 2 J2, J3 — > A _1 /iJ3, J4 — > A~ 2 /i _1 J4, J5 — > A _1 /i _2 J5, Jq — > \[i~ 1 Jq. 

There are four independent scaling-invariant monomial functions of the J„ , but it is convenient to consider 
the five functions given by u>x = J1J4, W2 — J2J5, W3 = J^Jq, W4 = J1J3J5, w$ = J2J4J6, which satisfy 
the single relation 

W1W2W3 = W4W5. (5.37) 

Since these Wi can also be expressed in terms of X\, . . . , x\q, the scaling- invariance implies that they can be 
written as functions of the symplectic coordinates yj as well. For instance, we have the formula 

(2/12/22/8 + 2/42/52/6 + 2/52/62/8 + 2/52/6)(2/i2/2 + 2/12/6 + 2/52/6 + 2/i + 2/6 ) 
Wl = , 

2/12/22/52/62/8 

and analogous formulae can be obtained for W2 — 0*wi, W3 — ((p*) 2 wi using the map (|5.36l) ; there are 
different expressions for W4 and w§ — ip*W4, but these are more unwieldy so they are omitted. Calculating 
the Poisson brackets between these functions, we have a subalgebra with four generators w±,W2, Ws,Wi, but 
this is more conveniently expressed with the extra function W5 included. The brackets can be determined 
from the relations 

{wi,w 2 } = wiw 2 -wi — w 5 , {wi, Wi} = wi(w 2 - w 3 ), (5.38) 

with all other brackets following by applying the map and noting that w\ cycles with period 3, and Wi has 
period 2, so (<^*) 3 wi = w\ and so (ip*) 2 Wi — Wi. The 5-dimensional algebra defined by (|5.38[) has three 
Casimirs, given by 

K.\ = 3 — wi — W2 — W3 + W5, K-2 = 3 — wi — W2 — W3 + Wi, C = W1W2W3 — 
where (as functions of J„) the quantities K,\ and /C2 come from (|5.22p . by taking the trace of the monodromy 

1/2 * A P 

matrix M„ for n — 1,2, while fixing C — corresponds to the constraint (|5.3T[) . The Casimir function 
B = JC\ + IC'2 = 3 — 27^i + H2, is also a first integral, with the components 

til = J\Ji + J 2 J5 + 'h-h = Wl +W 2 + W3, %2 = J1J3J5 + J3J4J6 = Wi + W 5 , 

which are homogeneous functions of the J„, are also first integrals, and Poisson commute: {^1,^2} = 0. 
Either of the latter two functions defines an integrable system on the two-dimensional symplectic leaves of 
the algebra generated by the Wi. The function C = ICyK.^ + B is another first integral that is also a Casimir. 
The functions K n , which cycle with period 4, are specified by the formula 

7- r Xji X n -|_g X n -\-i X n -\-g ^n+l^n+8 -, A 

K n = 1 1 1 1 , n = 1, ...,4. 

%n+6 Xn+3 2 ; n+ 2 2 ; n+6 Xn+3 x n+7 X n+2 X n +7 
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These quantities are invariant under the scaling (|5.30p . so they can also be written in terms of the iterates 
of Ij06]l : 

K _ Un+l + y n +6 + VnVn+l + Vn+lVn+6 + £fa+6j/n+7 _ -, . 

J\ n — , Ti — 1, . . . , 4, 

2M+32M+4 

where (p*y n — Un+i defines the sequence of y n for all n £ Z. Upon using (|5.35l) we find the relations 

{K U K 2 } = -K X K 2 + 1, {K U K 3 } = -K 2 + K 4 , 

which provide all the Poisson brackets between the K n by applying the cyclic property. Comparing with 
(|5.29[) and scaling by a factor of —3, this four-dimensional Poisson subalgebra is seen to be isomorphic to 
the algebra of the Ji in Example 15.161 Hence the Casimirs of this subalgebra are 

K\ = K 2 K 4 -K x - K 3 , &* 2 = KiK 3 - K 2 - K 4 , 

and, with B = —Hi + H 2 , two first integrals in involution are 

Hi = K 1 +K 2 + K 3 +K i , H 2 =KiK 3 + K 2 K 4 . 

From (|5.22l) . these Casimirs are shared with the subalgebra generated by the J n , and the two different 
formulae for B imply that the first integrals are related according to 

3-2H 1 +H 2 +H 1 -H 2 =0. 

Since {wi, Kj} = for all the functions B, C, Hi, Hi provide four first integrals that Poisson commute, 
as required for Liouville integrability of the 8-dimensional symplectic map (15.36[) . 



6 Integrable maps from Somos sequences 

The quadratic recurrences (|3.15l) (case (iv) of Theorem 13.121) , are referred to as three-term Gale- Robinson 
recurrences 16 , 34 . We mow consider the slightly more general case where these recurrences have coefficients: 

+ /3 (6.1) 

These can be included by adding extra nodes to the quiver for the coefficient-free recurrence (see Section 
10 in [13] )■ The coefficients are frozen variables, attached to these new nodes, which do not change under 
mutations. For three-term Gale-Robinson recurrences, one can add two extra nodes, corresponding to the 
parameters a, j3, so a quiver with N+2 nodes is obtained (Proposition 10.4. in [13]). It is then straightforward 
to check that the presymplectic form ui in (|2 . 10[) . defined by the same N x N skew-symmetric matrix (bjk) as 
for the coefficient-free case, is preserved by (|6.1[) . Hence Theorem 12.61 can be applied directly to the latter, 
to obtain a reduced symplectic map for suitable variables yj. 

Below we outline two different approaches to showing that the map for the variables yj is integrable. 
The first way is to use the fact that all of the recurrences (16.11) are ordinary difference equations that arise 
as reductions of the Hirota-Miwa equation, which is an integrable partial difference equation with three 
independent variables. The Lax pair of the Hirota-Miwa equation allows one to obtain Lax pairs for its 
reductions, and the associated spectral curves provide first integrals in terms of the yj. The second way is 
to find Somos- type recurrences of higher order that are satisfied by the iterates of (|6.1I) . The coefficients of 
these Somos-fc relations, for certain k > N, can also provide first integrals (analogous to the first integrals 
which appear as coefficients in linear recurrence relations for the iterates of the families (ii) and (iii)). 

Recently, Goncharov and Kenyon have found Somos recurrences arising as discrete symmetries of classical 
integrable systems associated with dimer models on a torus [19 . A further connection with relativistic 
analogues of the Toda lattice appeared in [7] . 
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6.1 Reductions of the Hirota-Miwa equation 

The relations (|6.ip all arise by reduction of the Hirota-Miwa (discrete KP) equation, which is the bilinear 
partial difference equation 

TilU = T 2 T_ 2 + T 3 T_ 3 . (6.2) 

In the above, T = T(ni, 712,7x3) is a function of three independent variables, and to denote shifts we have 
used T ±j = T\ rij ^ nj ±i. If we set 

T(ni,n2,na) = exp fa^SijUin^j r(n), (6.3) 

where S = (Sij) is a symmetric matrix and n = no + 5ini + 62TI2 + §3713, then x(n) satisfies the ordinary 
difference equation 

x(n + 6i)t(ti - 61) = ax(n + S 2 )x(n - 6 2 ) + (3x(n + 6 3 )x(n - 63) 

where a — exp(2(S , 22 — Sn)), j3 = exp(2(S33 — Sii)). Upon taking x n — x(n — 61) with 5i = ^N, 
S>2 = l(N- 2p), 63 = \{N -2q), this becomes (j6~T|) . 

In the combinatorics literature, equation (|6.2p is referred to as the octahedron recurrence, which has 
the Laurent property (shown in |llj). The Laurent property for three-term Somos (or Gale- Robinson) 
recurrences of the form (|6.1[) then follows by the reduction (|6.3|) . 

The Hirota-Miwa equation (|6.2p has a scalar Lax pair (see equation (3.8) in [24], for instance): it is the 
compatibility condition for the linear system given by 

r_ 1,3 ^1,2 + 7^2,3 = ^2,3^, / g4 N 

T -0-1,2 + T-1,3^2,-3 = r_i i2 ^, 

in terms of the scalar function -0 = ?/>(ni, 712, 713), with the same notation for shifts as before. Using the 
latter, one can use the reduction (16.31) to obtain Lax pairs for all of the Somos recurrences (|6.ip . which leads 
directly to spectral curves whose coefficients are conserved quantities. Here we briefly illustrate how this 
works for the cases N = 4 and N = 5 only; reducing the Lax pair becomes more involved as N increases. 
The general Somos-4 recurrence with coefficients is 

x n+i x n = ai„ +3 i„ + i + (ix 2 n+2 . (6.5) 

By taking the same monomials as in (|2.24p . this reduces to the map 

(2/1 , 2/2) ^ (2/2 , (ay 2 + P)/{yivl)) , (6.6) 

which becomes (|2.25p in the case a = /3 = 1, and preserves the same symplectic from uj. This means that 
(I6.6P has the invariant Poisson bracket 

{yi , 2/2} = 2/12/2- (6.7) 

The recurrence (|6.5p arises from the reduction (|6.3|) with 61 = 2, 62 = 1, 63 = 0. Upon taking Sjk = 
for j 5^ k, without loss of generality, and setting 



\ g •-'11 

2J rij log Aj + Sjjn 2 A x(n) <j>{n), with Ai = - -= , A 2 = AiC, A 3 = e 2S22 A 2 



the scalar linear equations (|6.4p reduce to a 2 x 2 linear system for the vector w = (<j>(n), <f>(n + I)) T . Up to 
shifts of indices, the coefficients of the latter can all be written in terms of the yj given in (|2.24p . as well as 
the spectral parameters ( and £, as follows. 
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Example 6.1 (Somos-4 Lax pair). The Lax pair for the map (|6.6p takes the form 

Lw = £w, w = Mw, (6-8) 

where the tilde denotes the shift n — > n+ 1. The matrices L = L(£), M = M(£) are functions of yj and the 
spectral parameter £, given by 



2/12/2 ' »/l2/2 1 M = I 

^C 2 -C (-yiift-^JC + iy' V-^kC^ 



2/11/2 ' 2/13/2 

The discrete Lax equation LM = ML holds if and only if the map (|6.6[) does. The spectral curve is 

det (L(C) - £ 1) = ( 2 + {H x ( - 1)£ + a 2 ( 3 + /3( 2 - 0, 
in which the coefficient of is the first integral 

H 1 =y 1 y 2 + — + — + (6.9) 
Vi 2/2 yiV2 

The level sets of Hi are biquadratic curves of genus one in the (2/1,2/2) plane, and the map is a particular 
instance of the QRT family [33 . 

Applying the reduction (|6.3[) to the Hirota-Miwa equation (|6.2p with S x = 5/2, 62 = 3/2, 5 3 = 1/2 leads 
to the general form of the Somos-5 recurrence with coefficients, which we denote by a, /3: 

(6.10) 

From the appropriate B matrix, one obtains j/i = x 1X4/(22 £3), 2/2 = 2:2 £5/ (23 £4) as coordinates in the plane 
(see [TS]), satisfying the Poisson bracket (|6.7p . and the corresponding map <p is also of QRT type, namely 



H> ■ (2/1 . 2/2) (2/2 , (52/2 + /3)/(2/i2/2))- (6.11) 



Similar calculations to those in the Somos-4 case yield the appropriate reduction of 

Example 6.2 (Somos-5 Lax pair). The map (|6.11l) arises as the compatibility condition LM = ML for a 
linear system of the form (16.8[) . where 

L = C + C 1 C + C 2 C 2 , M = C +f_^ °) C, (6-12) 

with Co=( !)• Cl =( 1 -uSt?ljL\ )' C2= l 5 + ^1 5 



The coefficient of in the equation for the spectral curve, that is 

det (L(C) - C 1) = e 2 - (25C 2 - JC + 1)C + 5 2 C 4 + /3C 3 = 0, 
gives a first integral whose level sets are cubic (also biquadratic) curves of genus one, that is 

J -2/1 +2/2 +5 f — + — ) (6.13) 
\Vi 2/2/ 2/12/2 
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Remark 6.3. The first integral (|6.9I) for Somos-4 can be rewritten in terms of the cluster variables, so 
that it becomes a ratio of homogeneous polynomials of total degree 4 in xi, X2, X3, £4, with the denominator 
just being X1X2X3X4. Similarly, in the case of Somos-5 the first integral (I6.13|) can be rewritten as a ratio 
of homogeneous polynomials of degree 5, with the denominator X1X2X3X4X5. It turns out that (|6.10|) has 
another rational first integral, also of degree 5 in terms of the cluster variables, which can be written as 

1 , 1 , 1\ , I ^ t. _ x i x m 



i = hhh + &[ T + T + T ) + -rrr with f 3 = ^^ for j = 1,2,3 

\Ji J2 hj J1/2/3 x j+1 

(see Proposition 2.3 in [3D])- However, the quantity / is not defined on the (2/1,2/2) plane, where there is only 
one first integral (as required for the Liouville- Arnold theorem). 

6.2 Bilinear relations of higher order 

In [32], Swart and van der Poorten proved that sequences generated by Somos-4 recurrences also satisfy 
quadratic (Somos-type) relations of order k, for all k > 4. They also noted that for Somos-5 sequences, there 

are Somos-fc relations of all odd orders k = 5, 7, 9, Moreover, the coefficients of the higher order relations 

are constant along orbits, which means that, as long as they are not trivially constant, they provide first 
integrals. 

In this subsection we explain how to obtain first integrals for Somos-7 recurrences using associated 
quadratic (bilinear) relations of higher order. The analogous results for Somos-6 recurrences are in [22] . 

To present the results concisely, it is convenient to consider the most general form of a Somos-7 recurrence, 
which is the four-term Gale-Robinson relation 

x n +7X n = (xx n+6 x n+ i + ftx n+5 x n+2 + y x n+i x n+3 . (6-14) 

With all three terms on the right hand side, this does not arise from a cluster algebra. Nevertheless, the 
general Somos-7 recurrence can be obtained as a reduction of the cube recurrence (Miwa's equation), and 
in [IT] this was used to prove the Laurent property for all four-term Gale-Robinson recurrences, including 
(|6.14p . In each of the cases where one of the parameters a, |3, y vanishes, (|6.14|) reduces to a bilinear relation 
with two terms on the right hand side, corresponding to a different cluster algebra in each case. 

Our main result on the family of recurrences (|6.14p is stated as follows and the remainder of this subsection 
is devoted to its proof. 

Theorem 6.4. The Somos-7 recurrence j6.14\ l has three independent first integrals, denoted Hi, T-L%, 1, which 
are rational junctions (in fact, Laurent polynomials) of degree 7 in Xi, . . . , £7. I can be written as 



I = P /1/2/3/4/5 + Y(/i/ 2 / 3 + / 2 / 3 / 4 + /aA/s) + ^{T 1 +T 2 +T 3 +T i +T, 

+ay (jhkh + ~hku + TUIh) + hhhhh + 137 \hhkfih + hfihhfs) + hfi fifth 



(6.15) 



in terms of the variables fj = XjXj + 2/x^ +1 , j = 1,2,3,4,5, while Hi and H2 are given in terms of yj — 
x i x j+'i/{ x j+iXj+2), as i n Proposition ESI below. For a = or y — 0, Somos-7 reduces to a Liouville 
integrable map tp in four dimensions, while for |3 = it reduces to an integrable map of the plane. 

The equation (|6.14p admits the action of the two-parameter scaling group (C*) 2 , via 

x n — > A /X™ x n (6.16) 

for non-zero A,/i. The variables fj are invariants under this scaling symmetry. In terms of these variables, 
the Somos-7 recurrence (|6 . 14[) is transformed to a recurrence of fifth order, namely 

fn+of n +Afn+'ifn+2fn+l fn = <X-fn+ifn+3fn+2fn+l + P fn+ifn+2 + J ■ (6-17) 
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Being of odd order, (|6.14l) has a further scaling symmetry depending on the parity of n: 

x n ->• v^ ir x n , veC*. (6.18) 

The variables f n have the symmetry /„ — > v ±2 f n for even/odd n respectively. The variables yj = fjfj+i arc 
the invariants under the combined scaling group and lead to 

yn+4Vn+3yl + 2yn+iy n = oiy n+3 y n+2 yn+i + (3 Vn+2 + Y- (6.19) 

Of the cluster algebra subcases, we shall see that the case a = and the case y = 0, the 4D map defined by 
(|6.19|) is symplectic, while for |3 = 0, additional scaling symmetries allow reduction to the plane. 

The quantities Hi,H2,I in Theorem 16.41 can all be written in terms of fj, so the recurrence (|6. 1 T[) has 
three independent first integrals. The rational function / is not invariant under (|6.18l) . which means that it 
does not provide a first integral for (I6.19|) . but both Hi and H2 do. 

The quantities Hi and H2 appear in coefficients of bilinear relations of higher order. Since the Somos-7 
recurrence is invariant under the three-parameter family of scalings defined by (|6.16p and (|6.18|) . one expects 
there to be relations of odd order having the same symmetry. The first non-trivial relation is the Somos-11 
recurrence (|6.20[) below. It can be seen that a combination of Hi and H2 appears in one of the coefficients. 
Since this coefficient remains constant along each orbit of (|6.14j) . it provides a non-trivial first integral. A 
second independent integral is provided by the Somos-13 recurrence (|6.2ip . given in the following. 



Proposition 6.5. The iterates of the Somos-7 recurrence j6.14\ l also satisfy the Somos-11 recurrence 

(6.20) 



Xn+u x n = - |3 x n+w x n+ i + y(y - (x 2 )x n+ $ x n+3 - a|3y x n+7 x n+i 
+ (a 5 + 2cty 2 + 2(3 3 + $H\ + (x 2 H 2 ) 



as well as the Somos-13 recurrence 

-|3y x n+ n x n+ 2 + Vx n+9 x n+4 + Ax n+8 x n+5 + Qx n+7 x n+6 , (6-21) 

where 

r = cc 4 y + cc 2 y 2 - a|3 3 + y 3 + otyH 2 , A = a 5 (3 + a 3 (3 y + 4a(3Y 2 + (3 4 + fi 2 Hi + (3 (a 2 + y)H 2 , 

8 = a 7 + 3a 3 Y 2 + 2a 2 (3 3 - ay 3 + a 2 |3"Hi + a 4 "H 2 
and, in terms of y.j = XjXj +3 / (x.j+iXj+2) for j = 1, 2, 3, 4, the first integrals Hi and H2 are given by 

Hi = Y2/12/22/32/4 + «P (2/12/3 + 2/12/4 + 2/22/4) + ay (yi + 2/2 + 2/3 + 2/4 + 
+a 2 (3 f^ + ^ + ^ + ^ + ^ + -*-) + (3y f- + - + - 



yiV3 i_ V2V4. 

V4 yi 



j_ 

V3 ' yi 



' \y1y2 V2ya J/3J/4 yiy3 V2yi J ' \y1y3 y2yi ) 1 \y\V2Vi y\V3Vi 

+a|3y ( + —^1— + ~Ar-) + ocy 2 ( — £ + , (6.22) 

1 1 \y1y2y3y4 yiyiyt yivivi ) ' \yiyiy3yi yiy2y£y4 J ' v ' 

n 2 = Y(2/i2/22/3 + 2/22/32/4) + a|3(2/i+ 2/2 +2/3 + 2/4) + ay (§ + f-) +a 2 |3 
+(3y f — + — + —) + y(a 2 + y) fr-V 

1 ' \y1y2 y2y3 V3Vi J lv '' \y1y2y3 



1V3 2/23/4 



- a p Y — + — — ) + «y; . 

1 ' \y1yiy3yi yiV2yiyi J yivzvivi 



Having obtained the first integrals, we now consider each of the three cases corresponding to cluster 
algebras separately, and explain how the reduction result of Theorem 12.61 applies in each case. 



The case a = 0: When a = 0, the recurrence (|6.14[) arises from a cluster algebra defined by a 7-node 
quiver. The latter comes from a 7 x 7 matrix, specified in terms of its columns by 

B = (-V3, -V4, vi + v 2 + 2v 3 + v 4 , -vi + v 4 , -vi - 2v 2 - V3 - v 4 , vi, v 2 ) , 
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where im B is spanned by the vectors 

Vj ■= e 3 ■- e j+1 - e j+2 + e j+3 for j = 1,2,3,4, (6.23) 
where is the jth standard basis vector. In this case, ker B is spanned by the three integer vectors 

u 1 = (1,1, 1,1, 1,1, 1) T , u 2 = (1,2,3,4,5,6,7) T , u 3 = (1, -1, 1, -1, 1, -1, if. (6.24) 
These three integer vectors produce a three-dimensional group of scaling transformations, 

x -> A U1 • ^ ■ ^ U3 • x, 

which coincides with the scalings defined by (16. 16)) and f|6 . 18[) . By Lemma 12.71 there is a symplectic form, 
given in terms of the scale-invariant monomials yj — x Vj for j — 1, 2, 3, 4, as 

dyi A dy 3 dy 2 A dy 3 dy 2 A dy 4 

LJ = 1 1 . 

2/12/3 2/22/3 2/22/4 

This yields the Poisson bracket 

{Vj> 2/j+i} = 0, {yj, 2/i+2> = -VjVj+2, {y j7 y J+3 } = VjUj+a- (6.25) 
In terms of one finds the map 

¥>■■ (2/1,2/2,2/3,2/4)^ (2/2,2/3,2/4, (02/3 +y)/(2/i2/22/32/4)), 

which is equivalent to iteration of (|6.19l) with a = 0, and preserves the nondegenerate Poisson bracket (|6 . 25|) . 
Setting a = in the two first integrals in (|6.22[) and computing their bracket gives {"Hi,7^2}|a=o = 0, so 
this is a Liouville integrable system in 4D. 

The case |3 = 0: In this case the matrix B is specified by B — (— V2, Vi,V2, — vi + V2, — Vi, — V2, Vi), where 
Vj = Vj + Vj+i + v J+ 2 for j — 1, 2. The kernel of B is 5-dimensional, being spanned by the integer vectors 
u 2 and u 3 together with u 4 = (1,0, 0, 1, 0, 0, 1) T , u 5 = (0, 1, 0, 0, 1, 0, 0) T , u 6 = (0, 0, 1, 0, 0, 1, 0) T . These 
five independent integer vectors give an action of the algebraic torus (C*) 5 on x by scaling transformations, 
and the scalings (|6.16|) and (|6.18| form a three-parameter subgroup, since ui = 114 + 115 + u@. The invariants 
under the full 5-parameter scaling group are 

Vj = x Vj = XjXj+s/^Xj^Xj+a) = yjyj+iyj+2, 

and in terms of these one obtains the general Lyncss-2 recurrence with coefficients, that is 

2M+2 y n = oc y n +i + y, (6.26) 

which is equivalent to iteration of a map in the (j/i, 7/2) plane and of QRT type [33j . with invariant symplectic 
form (2/1 i^fMj/i A d{/2- The first integral Hi does not reduce to the plane, as it is not invariant under the 
full 5-dimensional scaling group, but H2 is fully invariant, and reduces to a first integral of (|6.26[) : 

n 2 \^o = y (2/1 + 2/2) + ay ( |i + ^ + y(a 2 + y) ( I + 1) + ^ 



,2/2 2/1/ V2/1 2/2 J 2/12/2 

The level sets of the latter are cubic (and biquadratic) plane curves of genus one. 

The case y = 0: This is very similar to the case a = 0, since the recurrence comes from the rank 4 matrix 

B = (— V2 - V4, Vl + V 2 + V4, -vi + V 2 , -V2 + V 3 , -V3 + V 4 , -Vl — V3 — V4, Vl + v 3 ) , 
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with the same vectors Vj as in (|6 . 23|) . Since both ker B and imB are the same as for a = 0, there is the 
same scaling group with invariants yj = x Vj . The symplectic form in this case is 

dyi A dj/2 . dyx A dy 4 . dy 3 A dy 4 

u> = 1 1 , 

V1V2 2/12/4 2/32/4 

and (up to an overall constant) this gives the unique log-canonical Poisson bracket 

{vj,yj+i} = VjVj+u {vj,Vj+2} = o = {yj,yj+a} (6.27) 

that is preserved by the map 

<P- (2/1,2/2,2/3,2/4)^ (2/2,2/3,2/4, (at/22/4 + P)/(2/l2/22/32/4))- 

The latter map in 4D corresponds to iteration of (|6. 19|) with y = 0. Setting y = in (|6.22|1 and computing 
the bracket using (I6.27|) gives {"Hi, %2}|y=o = 0, so the two first integrals are involution, as required. 

Remark 6.6. We have shown that for the two parameter subcases of (|6.14[) . the reduced map is integrable 
in the Liouville sense, either in 4D or 2D. However, there is a different Poisson structure in each case. 
When oc|3y 7^ 0, it is easy to check that there is no log-canonical Poisson bracket in the variables yj that is 
compatible with f)6. 19[) . Nevertheless, we expect that there is a compatible Poisson structure for which the 
two first integrals in (I6.22j) are in involution, so that (|6.19[) defines a 4D map that is Liouville integrable. 
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